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1.  INTRODUCTION 


The  on-going  effort  whose  various  facets  have  been  described 
e.g.  in  [Blaha,  1975,  1977,  1979,  1980,  1981]  is  concerned  with  the  adjust¬ 
ment  of  satellite  altimeter  data  in  two  steps,  the  first  performed  in  terms 
of  a  truncated  set  of  spherical-harmonic  (S.H.)  potential  coefficients  and 
the  second  performed  in  terms  of  point-mass  (P.M. )  magnitudes  as  parameters. 
The  emphasis  has  been  shifted  recently  toward  SEASAT  altimeter  data  as  the 
main  source  of  observed  quantities  in  the  first  adjustment;  gravity  anoma¬ 
lies  and  other  sources  of  geopotential  information  have  been  included  via 
the  weighted  S.H.  coefficients.  At  this  stage  (first  adjustment),  six 
weighted  state  vector  (s.v.)  parameters  per  orbital  arc  are  also  included 
in  the  simultaneous  least-squares  process.  The  first  adjustment  is  global 
in  character.  One  of  its  most  important  products  is  a  revised  set  of  S.H. 
coefficients  which  may  be  of  interest  in  itself,  and  which  is  especially 
useful  in  predicting  geoid  undulations,  gravity  anomalies  and  other  quanti¬ 
ties  related  to  the  disturbing  potential  (such  as  deflections  of  the  vertical 
or  gravity  gradients)  on  the  global  scale. 

The  data  for  the  second  adjustment  consist  mainly  of  the  residuals 
from  the  first  adjustment,  although  other  quantities  (gravity  anomalies, 
deflections  of  the  vertical,  etc.)  can  enter  this  phase  independently.  In 
this  process  a  more  detailed,  but  regional,  geoid  is  derived.  Predictions 
of  the  other  quantities  just  mentioned  can  also  be  made  in  the  region  of 
interest.  A  given  set  of  point  masses  has  a  chosen  distribution  which,  as 
a  rule,  is  uniform  and  is  characterized  by  the  1.6:1  depth-side  ratio. 
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Denser  sets  of  point  masses  could  be  superimposed  on  the  basic  set,  leading 
to  an  even  more  detailed  description  of  the  gravity  field  in  specific  sub- 
regions.  This  can  be  related  to  the  mean  values  of  geoid  undulations 
derived,  for  blocks  of  a  certain  size,  from  satellite  altimetry.  In  some 
areas  of  pronounced  and  varied  geoidal  relief,  such  as  in  the  Puerto  Rico 
trench  area,  the  differences  between  the  mean  and  the  actual  undulations 
could  become  large  and,  thus,  smaller  blocks  might  be  chosen  to  describe 
the  geoid.  From  geoid  undulations  (or  their  means)  one  could  derive  other 
quantities  related  to  the  disturbing  potential,  upon  using  the  appropriate 
cross-covariance  functions.  The  present  approach  with  point  masses  cir¬ 
cumvents,  by  construction,  the  need  for  these  functions. 

One  final  product  obtained  with  the  P.M.  adjustment  superimposed 
on  the  S.H.  adjustment  is  a  set  of  contour  maps  based  on  predicted  values 
for  geoid  undulations,  gravity  anomalies,  etc.  A  typical  example  is  a 
geoidal  map  in  a  region  containing  point  masses.  This  region  exhibits  detail 
ed  geoidal  features,  while  far  from  it  the  geoid  described  by  the  potential 
coefficients  alone  is  very  smooth;  the  transition  from  one  region  to  the 
other  is  gradual. 

Since  the  global  S.H.  adjustment  is  the  basis  for  any  detailed 
geoidal  resolution,  it  has  been  refined  on  several  occasions,  both  from  the 
accuracy  and  the  economy  standpoints.  For  example,  criteria  have  been  es¬ 
tablished  specifying  the  maximum  and  the  minimum  allowable  length  of  satel¬ 
lite  arcs.  In  terms  of  SEASAT  altimetry  these  criteria  translate  into  seven 
minutes  and  about  one  minute  in  duration,  respectively,  as  described  in 
Chapter  2  of  [Blaha,  1981].  The  same  report  contains  also  the  development 
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resulting  in  a  substantial  improvement  in  the  economy  of  the  short-arc  j 

algorithm,  achieved  through  a  reduction  in  the  number  of  S.H.  potential  j 

coefficients  entering  the  orbital  integrator.  The  corresponding  reduction  \ 

in  the  run-time  requirements  does  not  compromise  the  quality  of  the  alti-  j 

metric  information.  Indeed,  the  information  represented  by  the  SEASAT  j 

altimeter  data,  the  NSWC  precise  ephemeris,  the  GEM  10  S.H.  potential  coef-  j 

ficients  (at  least  within  the  16,16  truncated  set),  as  well  as  the  reference  j 

field  parameters  has  been  found  to  be  of  excellent  quality.  The  pertinent  j 

analysis  is  presented  in  Appendix  1  which  draws  a  tentative  conclusion  with 
regard  to  the  closed- form  expression  giving  the  degree  variances  for  geoid 
undulations. 

Due  to  the  high  quality  of  the  SEASAT  altimetric  system  compared  I 

■j 

to  its  predecessors  (e.g.,  the  GE0S-3  system  including  the  broadcast  ephemeris), 
modeling  errors  caused  by  certain  sea  level  changes  can  no  longer  be  ignored.  j 

The  most  important  changes  are  those  caused  by  the  tide-generating  forces  of 
the  moon  and  the  sun.  Of  these,  eleven  long-period,  diurnal  and  semidiurnal 
tidal  effects  are  discussed  in  the  body  of  this  report,  starting  with  the 
equilibrium  tide.  The  initial  development  overlaps,  to  a  certain  extent, 
with  Sections  3.1  and  3.2  of  [Blaha,  198l]  where,  however,  fewer  tidal  con¬ 
stituents  were  treated.  ( 

In  addition  to  this  Final  Report,  the  reports  covering  the  entire 
period  of  the  present  research  contract  have  been  [Blaha  1979,  1980,  1981].  j 
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2.  EQUILIBRIUM  TIDE 


When  determining  the  height  of  the  theoretical,  or  equilibrium, 
tide,  its  individual  component  "h."  is  associated  with  the  tidal  con- 

J 

stituent  "j"  of  amplitude  A.  and  argument  a..  The  total  height  is  then 

J  J 

the  sum  of  the  individual  h ^ ' s .  The  basic  formulas  adopted  in  this  de¬ 
velopment  are  (129)  -  (131)  of  [USCGS,1958]  abbreviated  here  as  [US] . 

Since  only  the  average  values  of  the  h.'s  with  regard  to  the  longitude 

J 

of  the  moon's  mode  are  sought  at  a  first  stage,  the  "node  factor",  f,  is 
taken  as  unity  in  the  pertinent  formulas.  This  implies,  for  example,  that 
the  "permanent  tide"  symbolized  by  Aq  is  considered  through  its  mean  effect 
over  one  full  revolution  (or  several  full  revolutions)  of  the  moon's  node; 
such  a  revolution  is  completed  in  about  18.6  years. 

In  considering  the  average  (in  the  above  sense)  combined  effect 
of  the  moon  and  the  sun,  the  constituent  height  can  be  expressed  by 

h.  =  A.  cos  a.  .  (2.1) 

J  J  J 

The  amplitude  varies  with  4>,  the  geocentric  latitude,  as 

Aj  >  Kj  Gj(»)  Cj  ;  (2.2) 

here  f  =  l  is  assumed  so  that  the  coefficient  C.  represents  the  mean  value 

J 

of  a  Dertinent  function  with  respect  to  the  longitude  of  the  moon's  node. 

In  the  following,  three  indices  (a,  b,  c)  will  be  used: 

a. ..  long-period  constituents, 

b. ..  diurnal  constituents, 

c. ..  semidiurnal  constituents. 
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The  meaning  of  the  symbols  K.  and  G.(<j>) 

J  J 

is  thus  narrowed  down  to 

K,  =  %Ga  =  0. 13334  m  , 

(2.3a) 

K.  =  K  =  Ga  =  0.2667  m  , 
b  c 

(2.3b) 

where 

G  =  (3/4)(M/E)(a/rz)  , 

(2.3c) 

with  M  and  E  being  the  moon's  and  the  earth's  masses,  respectively,  a  being 

the  mean  radius  of  the  earth  (6,371km)  and  rM  being  the  mean  earth-moon 
distance;  according  to  the  values  listed  in  Table  1  of  [US] ,  G  -  0.41865  x  10 

Further, 


Ga(<J>)  =  (1-3  si n24>)  ,  (2.4a) 

Gb(<t>)  =  sin  2<J>  ,  (2.4b) 

G  (<f>)  =  cosz4>  •  (2.4c) 

v 

In  [Blaha,  1981],  Table  1  was  constructed  featuring  the  constitu¬ 
ent  heights  expressed  according  to  (2.1)  -  (2.4c).  Its  "extreme  magnitude" 
column  features  the  largest  values  of  h.  which  can  be  reached  as  follows: 

J 

a  ...  <f>  =  ±90°  , 
b  ...  $  =  ±45°  , 
c  . . .  =  0  . 

In  each  group  (a,  b,  or  c)  the  constituents  are  listed  in  the  descending 
order  of  magnitude.  This  table  is  reproduced  here  for  the  sake  of  complete¬ 


ness. 
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SYMBOL 


DESCRIPTION 


EXTREME 

MAGNITUDE 


Ao 

Mf 

Mm 

SSa 

constant 

semimonthly 

monthly 

semiannual 

.7384 

.1566 

.0827 

.0728 

.0985m  (1-3  sin2<f>) 

.0209m  (1-3  sinz<j))cosaMf 
.0110m  (1-3  sin2<j>)cosa^m 
.0097m  (1-3  sin2<|>)cosaSSa 

.197m 

.042m 

.022m 

,019m 

K1 

declinational 

1  uni -solar 

.5305 

.1415m  sin2<j>  cosai/ 

*1 

.  141m 

°1 

principal  lunar 

.3771 

,1006m  sin2cj>  cosotQ^ 

.101m 

P1 

principal  solar 

.1755 

.0468m  sin24>  cosap^ 

.047m 

m2 

principal  lunar 

.9085 

.2423m  cos2cj>  cosa^ 

.242m 

S2 

principal  solar 

.4227 

.1127m  cos2c|>  cosac 
^2 

.  113m 

n2 

ecliptical  lunar 

.1759 

.0469m  cos2<{>  cosa^ 

.047m 

k2 

declinational 

luni-solar 

.1151 

.0307m  cos2<J>  cosou, 
k2 

.031m 

Table  1 

Approximate  heights  of  selected  tidal  constituents, 
including  extreme  magnitudes 


2.1  Long-Period  Constituents 

If  only  the  average  effects  were  of  interest,  the  constituent 
heights  h^  in  the  group  "a"  of  Table  1  could  be  compared  with  the  formulas 
of  [Lisitzin,  1974],  pages  49  and  38,  39,  respectively,  implying  for  the 
tidal  potential  of  the  long-period  terms: 

Wfl  =  0.96621  (1-3  sin24>)  m2/sec2, 

o 

WM^  =  0.20460  (1-3  sin2<|>)  cos  2s  m2/sec2, 

=  0.10796  (1-3  s i n 2 )  cos(s-p)  m2/sec2, 

WSSa  =  °*09531  (1_3  sin 2<f>)  cos  2h  mVsec2, 

where 

h  =  mean  longitude  of  the  sun, 
s  =  mean  longitude  of  the  moon, 
p  =  longitude  of  the  lunar  perigee. 

It  follows  for  the  constituent  heights: 

Wfl  /  g  -  0.0986m  (1-3  sin2ij>)  , 

Ao 

wMf/ ?  -  0.0209m  (1-3  sin2<j>)  cos  2s  , 

WMn/  9  ~  O-OllOm  (1-3  sin24>)  cos(s-p)  , 

Ws$a/g  ~  0. 0097m  (1-3  sin2^)  cos  2h  , 

where  "g,  the  average  terrestrial  gravity,  is  adopted  as  9.80  m/sec2.  Except 
for  the  less  explicit  notations  for  the  arguments  "a"  in  Table  1,  the  agree¬ 
ment  is  nearly  perfect. 
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If  a  more  accurate  representation  is  sought,  the  influence  of 
the  moon  and  the  sun  must  be  treated  separately.  For  this  purpose,  the 
node  factor  for  the  moon,  denoted  in  general  as  f.,  can  no  longer  be  assumed 

J 

to  be  unity  (it  is  always  unity  for  the  sun).  When  applied  separately  for 
the  moon  and  the  sun,  (2.2)  becomes 


ft3  ’ 

K,  G.(4>)  C.  f.  ...  moon  , 

J  J  J  J 

(2.5a) 

a;  ■ 

K.  G.(4>)  Cl  ...  sun  . 

J  J  «J 

(2.5b) 

The  node  factor  (a  function  of  the  longitude  of  the  moon's  node  with  the 

periodicity  of  about  18.6  years)  changes  very  slowly  from  year  to  year  for 

each  constituent.  Table  14  of  [US]  gives  the  value  of  the  pertinent  f. 

3 

for  the  middle  of  each  year  between  1850  and  1999. 

If  a  given  constituent  represents  the  moon's  action  alone,  (2.1) 
is  adopted  without  change  in  notations  and  A.  is  computed  according  to 

J 

(2.5a).  For  a  strictly  solar  constituent,  A.  and  a.  in  (2.1)  are  replaced 

\) 

by  A',  and  a'.,  with  A',  computed  as  in  (2.5b).  If  a  constituent  is  composed 
3  3  3 

of  both  effects,  the  resulting  h.  is  obtained  as 

(2.6) 

h.  =  A.  cosa.  +  Al  ccsal  =  K .G . (<J>) (C .f .  cosa.  +  Cl  cosa'.). 

33  33  333  33  33  3 


For  the  permanent  tide,  the  value  which  corresponds  to  the  active 
life-span  of  SEASAT  can  be  associated  with  mid-year  1978  and  is  given  as 


Mm 


1.131  . 


(2.7) 
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Since  a.  and  a',  are  immaterial  for  the  permanent  tide,  with  the  aid  of 

J  J 

(2.3a)  and  (2.4a)  equation  (2.6)  becomes 

hA  =  0.13335m  (1-3  sin2<J,)(CA  fft  +  C ^  )  .  (2.8) 


Table  2  of  [US]  gives 


CA  =  0.5044, 


«  0.2340  , 


which,  when  added  algebraically,  yield  the  value  0.7384  seen  in  Table  1 

for  an  average  effect.  However,  in  considering  the  specific  case  of 

SEASAT  altimetry  and  the  corresponding  value  f.  in  (2.7),  equation  (2.8) 

Ao 

yields  the  explicit  form  for  the  equilibrium  height  of  the  constituent  A  : 


hfl  =  0.1073m  (1-3  sin24>)  . 


(2.9) 


The  constituents  Mf  and  Mm  are  due  exclusively  to  the  moon.  In 
agreement  with  [US]  the  tidal  arguments  are  written  as 


aMf  =  aMf  "  2?  ’ 

“Mm  =  “Mm  ’ 

Table  6  or  Table  11  of  [US]  give  for  the  SEASAT  observational  epoch: 
-2C  =  -1.07°  . 

The  coefficients  "C"  and  the  node  factors  "f"  for  this  epoch  are 

=  0.1566  ,  f|^  =  0.625  , 

CMm=  0-0827  ’  fMm=  L131  * 
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According  to  (2.1)  and  (2.5a),  the  equilibrium  formulas  representing  the 
SEASAT  epoch  thus  read 


hMf  =  0.0131m  (1-3  si n24>)  cos  (c^f  -  1.07°)  , 
hMm  =  0.0125m  (1-3  sin24>)  cosa^  . 

The  constituent  SSa  owes  its  existence  to  the  sun  only, 
and  (2.5b)  apply.  According  to  [us] ,  one  can  write 

aSSa  =  aSSa  ’ 

CSSa  =  °'0728  » 

yielding,  regardless  of  the  SEASAT  epoch, 

hSSa  =  °-00^7m  C 1-3  sin2d>)  cosa£Sa  , 


(2.10) 

(2.11) 

hence  (2.1) 


(2.12) 


which  is  essentially  (except  for  the  argument  notation)  the  same  as  the  cor¬ 
responding  expression  in  Table  1. 


2.2  Diurnal  Constituents 


When  considering  the  constituent  K1  ,  one  realizes  that  the  value 

h„  is  made  up  of  the  moon's  and  the  sun's  contribution,  hence  it  is  of  the 
Ki 

form  (2.6).  Due  to  the  changing  longitude  (N)  of  the  moon's  node,  a„  differs 

K1 

from  a!,  by  a  small  quantity  -v,  namely 
K1 

a„  =  a!,  -  v  .  (2.13) 

h  h 

Similar  to  the  node  factor,  the  periodicity  of  v  is  approximately 
18.6  years.  Over  short  periods  (e.g.  less  than  a  year)  it  can  be  considered 
constant.  Table  6  of  [US]  gives  the  values  of  v  according  to  N.  For  the 
beginning  of  September  1978,  the  epoch  which  is  quite  representative  of  the 
SEASAT  data  series,  the  value  of  N  found  from  Table  4  of  [US]  is 

N  =  178°  ,  (2.14) 


implying  that 

v  «  0.57°  .  (2.15) 

The  value  in  (2.14)  is  considered  constant  for  all  SEASAT  observations. 

In  order  to  obtain  h„  from  (2.6)  in  the  form  similar  to  (2.1)  in 

K1 

conjunction  with  (2.5a),  one  has  to  simplify  the  following  expression: 

cv  =  Cv  fv  cos(a^  -v)  +  Ci  cosot|J,  ,  (2.16a) 

which  corresponds  to  the  quantity  inside  the  parentheses  in  (2.6)  with  Kj 
substituted  for  j  and  with  (2.13)  taken  into  account.  Equation  (2.16a)  can 


jKSjk 

i 

i*  * 

be  developed  into 

:  i 

V 

(C„  f„  cosv  +  Cw  )  cosa!,  +  Cw  fw  sinv  sinaw 

,  (2.16b) 

'  i 

which  is  of  the 

form 

•i 

c  =  c1 

cosa  +  C2  sina  =  (C^  +  [  (C^  cosa  +  C2  sina)/(C^  +  C^) “]. 

» 

i 

*  *! 

If  v’  is  defined 

as  arc  tg  (Cg/C^),  one  has 

k.  1 

sinv' 

=  C2/(C\  +  q)h  ,  cosv1  =  c1/(c|  +  C|)H  , 

and  hence 

t 

C  = 

(Cj  +  C\Y2  cos  (a  -  v* )  . 

i 

; 

When  applied  to 

(2. 16b)  this  yields 

i 

CK 

Ki 

(C„  fv  )  cosa,,  , 

K1  K1  h 

(2.17a) 

where 

i 

1 

1 

VK1 

\ 

=  (Cj  +  cl)55. 

(2.17b) 

i 

i 

< 

\ 

1 

C1  = 

CK  fK  cosv  +  ,  C2  =  CK  sinv  ’ 

(2.17c) 

a.,  — 

h 

otjJ,  -  v1  ,  v'  =  arc  tg  ^/Cj)  . 

(2. 17d) 

i 

.  t 

The  coefficient  Cv  represents  a  certain  mean  value 

h 

defined  as 

> 

~CK  = 
K1 

mean  [  (Cj  +  C^)*5  cosv']  =  mean  , 

(2.18a) 

i 

' 

► 
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U.  fc 

i 

which  in  Table  2  of  [US]  is  listed  to  be 


L  =  0.5305  .  (2.18b) 

*1 

Further  listed  are  the  values 

C„  =  0.3623  ,  C*  =  0.1681  . 

K1  K1 

The  node  factor  for  K1  is  identical  to  that  for  Jj  of  equation  (76)  of 
[US],  and  it  is  listed  for  mid-1978  as 

f„  =  f,  =  0.827. 

K1  J1 

With  these  values  for  C„  ,  CA  ,  f„  ,  and  with  v  from  (2.15),  one  calculates 

i  K1  K1 

(see  2.17  c,d): 

Cl  *  0.4677,  C2  »  0.002981,  v'  =  0.37  ; 
from  here  it  follows  (see  2.17  b,d)  that 

C„  fy  =  0.4677  ,  cL  «  a'  -  0.37°.  (2.19) 

K1  K1  h  K1 

The  values  in  (2.19)  could  also  be  obtained  more  directly  from 
Table  14,  and  Table  6  or  Table  11  of  [US] ,  respectively.  In  particular, 
0.4677  could  be  found  upon  multiplying  0.5305  in  (2.18b)  by  the  corresponding 
node  factor  ?K  ,  listed  for  mid-1978  in  Table  14  (under  the  heading  K^)  as 
0.882;  and  0.37°  could  be  found,  for  N=178°,  from  Table  6  (under  the  heading 
v')  or  from  Table  11  (under  the  heading  K^).  In  either  case  the  equilibrium 
formula  (2.6),  applied  to  Kj  in  conjunction  with  (2.3b),  (2.4b),  (2.17a)  and 
(2.19),  becomes 

h„  =  0. 1248m  sin2<J>  cos(ai  -0.37°)  , 


(2.20) 


referring  to  the  epoch  of  SEASAT  altimeter  data  acquisition.  The  variable 
part  of  the  argument  for  this  and  other  tidal  constituents,  a^,  will  be 

«J 

described  later. 

The  effect  of  0^  is  due  exclusively  to  the  moon,  hence  (2.1) 

and  (2.5a)  apply.  In  analogy  to  (2.13)  and  the  development  that  followed, 

an  will  be  written  as  al  (this,  in  itself,  is  immaterial  here)  plus  some 
U1  U1 

small  quantity  which  will  again  be  considered  constant  due  to  the  short 
active  life-span  of  SEASAT.  In  particular, 

cu  =  cd  +  (2? -v)  .  (2.21) 

^1  U1 

For  N  given  in  (2.14),  Table  4  of  [US]  yields  approximately  0.54°  for  £ 
and  0.57°  for  v  (see  equation  2.15  above);  Table  11  of  the  same  reference 
yields  directly  2£  -  v  under  the  heading  0^  In  either  case  the  result  is 

2£  -v  *  0.50°.  (2.22) 

The  coefficient  "C"  and  the  node  factor  for  mid- 1978  are 

Cn  =  0.3771,  fn  =  0.806.  (2.23) 

U1  U1 

Upon  inserting  the  results  (2.3b),  (2.4b)  and  (2.21)  -  (2.23)  into  (2.1)  and 
(2.5a),  one  has  the  equilibrium  formula  for  0^  representing  the  SEASAT  obser¬ 
vational  epoch: 

hg  =  0.0811m  sin2<t>  cos(oig  +  0.50°)  .  (2.24) 
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extract 


The  constituent  Pj  is  due  to  the  sun's  effect.  From  [US]  we 

ap  =  ap  » 

P1  K1 

C0  =  0.1755  , 

K1 

yielding 

hD  =  0.0468m  sin24>  coscxi  .  (2.25) 

K1  1 

This  expression  is  independent  of  the  SEASAT  epoch  and  is  essentially  identi¬ 
cal  to  the  corresponding  expression  in  Table  1. 


2.3  Semidiurnal  Constituents 


The  most  important  of  all  the  tidal  constituents.  M2>  is  due 
exclusively  to  the  moon.  It  can  be  developed  in  a  complete  analogy  to  the 
approach  followed  for  0r  The  argument  is  written  as 

aM  =  aM  +  (2£-2v)  ,  (2.26) 

2  2 

where  £  and  v  were  already  found;  thus 

25  -  2v  =  -0.07°  ,  (2-27> 

which  is  also  the  value  given  in  Table  11  of  [US]  under  the  heading  M?. 
Further,  one  has 


CM  =  0.9085  , 
m2 

fM  =  1.038  , 
m2 

(2.28) 

hu  =  0.2515m  cos2<t> 

cos(a^  -  0.07°)  , 

(2.29) 

which  is  the  equilibrium  formula  for  corresponding  to  the  SEASAT  obser¬ 
vational  epoch.  It  has  been  obtained  from  (2.26)  -  (2.28)  in  the  same  way  as 
(2.24)  was  obtained  from  (2.21)-  (2.23)  except,  of  course,  that  (2.4b)  has 
been  replaced  by  (2.4c). 

The  constituent  S2  owes  its  existence  to  the  sun.  The  argument 

is  thus  written  as  aA  in  agreement  with  the  original  convention,  and  the 

b2 

node  factor  is  omitted.  In  other  respects  the  equilibrium  formula  for  S2 
is  derived  similar  to  (2.29)  above,  namely 

hc  =  0.1127m  cos 2 <p  cosal  .  (2.30) 

S2  b2 
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In  this  case,  no  special  considerations  related  to  the  SEASAT  observational 
epoch  are  necessary. 


The  constituent  N2  is  due  solely  to  the  moon;  one  has 

aN2  =  +  {2?  ‘  2v)  ’ 

where,  according  to  (2.27), 

2?  -  2v  =  -0.07°  . 


Furthermore,  one  finds 


hence 


CN  =  0.1759, 
n2 


fN  =  1-038  , 
in2 


hN  =  0.0487m  cos2<|>  cos(c^  -  0.07°)  ,  (2.31) 

which  is  the  equilibrium  formula  for  N2  corresponding  to  the  SEASAT  epoch. 

The  constituent  K2  reflects  the  effect  of  both  the  moon  and  the 
sun,  hence  (2.6)  applies,  where 


Accordingly, 


2v  . 


h„  *  Ga  cos2<J>  c„  , 

*2  i'2 

where  Ga  was  given  in  (2.3b)  and  where 


C  f 

lk2tk2 


cos(<xJ,  -  2v)  +  C!.  cos  aj, 
*2  N2  *2 
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1 


in  analogy  to 

(2.16a).  The 

numerical 

values  are 

2v 

=  1.14°  ; 

CK 

c 

=  0.0786  , 

fK  = 
c 

0.630  ; 

CK 

K2 

=  0.0365  . 

From  this  point  on,  the  algorithm  contained  in  the  formulas  (2- 16a)- (2. 19) 


can  be  adopted  in  its  entirety,  except  for  the  following  changes  in  notations: 
l<i  Kg  ,  v  -*■  2v  ,  v1  -  2v"  . 


The  results  are 


(C„  f„  )  cosoL 
*'2  ■'2 


L  f„  =  0.08602  , 
1^2  ^2 


-  2v"  , 


2v"  =  0.66°  . 


Thus 


h.,  =  0.0230m  cos24>  cos(a.)  -  0.66°)  . 

^2  In  g 


(2.32) 


This  outcome  could  also  be  obtained  from  Tables  2,  14,  and  6  or  11  of  [US]. 

In  particular,  Table  2  lists  (L  (see  Note  4  on  page  166)  as  0.1151.  The 

.  K2 

corresponding  node  factor  f„  is  0.748  appearing,  for  the  mid-1978,  in  Table 

k2 

14;  this  yields  the  product  0.08609  which  is  in  a  good  agreement  with  the 
above  computed  value.  Table  6  gives  2v"  for  N =  178°  (see  equation  2.14)  as 
0.66“;  equivalently.  Table  11  shows  this  value  under  the  heading  K^. 
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2.4  Tidal  Arguments 

The  explicit  expressions  for  the  equilibrium  tidal  arguments 
are  developed  in  a  way  similar  to  Table  2  of  [US]  with  a  few  minor  changes. 

One  change  pertains  to  the  constant  part  "u"  which  is  presently  expressed 
numerically  (in  °)  and  represents  the  SEASAT  observational  epoch.  With 
regard  to  the  computation  of  at,  the  variable  part  of  the  argument  (in 

J 

Table  2  of  [US]  denoted  as  V),  UT  ±12  hours  is  used  instead  of  T,  the  hour 
angle  of  the  mean  sun  at  Greenwich  at  the  time  of  the  tidal  evaluation. 

Since  all  the  quantities  will  be  considered  as  given  in  degrees  instead  of 
hours,  the  following  applies: 

T  =  UT  ±180°  ,  (2.33) 

where  UT  (in  0 '  is  obtained  by  multiplying  UT  (in  hours)  by  15  (°/hour), 
etc.  The  other  two  variables  needed  for  the  evaluation  of  a",  at  Greenwich 
for  the  desired  constituents  are  h,  the  mean  longitude  of  the  sun,  s,  the 
mean  longitude  of  the  moon  and  p,  the  longitude  of  the  lunar  perigee.  In 
terms  of  local  —  rather  than  Greenwich  —  arguments,  UT  is  replaced  by 
UT  +  X,  where  X  (in  °)  is  the  customary  east  longitude  of  the  point  where  the 
tidal  evaluation  is  sought,  symbolized  by 

local  argument  ...  UT-*-  UT  +  X  .  (2.34) 

In  order  to  indicate  the  computation  of  the  equilibrium  tidal  arguments  at 
Greenwich,  Table  2  lists  these  arguments  in  two  parts  (see  its  second  column), 
a",  and  u.;  the  final  argument  is 

J  J 

a.  =  at  +  u.  .  (2.35) 

J  J  J 
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CONST I- 

GREENWICH  ARGUMENT 

aj 

PERIOD 

TUENT 

“j 

+u  • 

J 

o/day 

Jan  0.5,  1900 

Mf 

2s 

-1.07° 

26.352793539 

1.09803306 

180.874844° 

13.6608  d 

Phi 

s  -  p 

+0.00 

13.064992739 

0.54437470 

296.109403° 

27.5546  d 

SSa 

2h 

+0.00 

1.971294671 

0.08213728 

199.393356° 

182.6211  d 

K1 

UT  +  h  +  90° 

-0.37° 

360.985647335 

15.04106864 

189.696678° 

23.9345  hr 

°1 

UT  -  2s  +  h  -  90° 

+0.50° 

334.632853797 

13.94303557 

188.821833° 

25.8193  hr 

P1 

UT  -  h  -  90° 

+0.00 

359.014352665 

14.95893136 

170.303322° 

24.0659  hr 

_ 

m2 

2UT  -  2s  +  2h 

-0.07u 

695.618501132 

28.98410421 

18.518511° 

12.4206  hr 

S2 

2UT 

+0.00 

720. 

30. 

0. 

12  hr 

n2 

2UT  -  3s  +  2h  +  p 

-0.07° 

682.553508393 

28.43972952 

82.409108° 

12.6583  hr 

k2 

2UT  +  2h 

-0.66° 

721.971294671 

30.08213728 

199.393356“ 

11.9672  hr 

Table  2 

Greenwich  arguments  and  related  quantities 
for  selected  equilibrium  tidal  constituents 


The  aj  part  agrees  with  [Schwiderski ,  1980],  page  172,  and  with  [Lisitzin, 
1974],  page  12. 

For  the  explicit  computation  of  al,  the  expressions  for  h,  s  and  p 

J 

are  adapted  from  [US],  page  162,  as 
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h  =  279.696678°  +  36 ,000. 768925° T  +  0.000303°T2  ,  (2.36a) 

s  =  270.437422°  +  48 1,26 7. 892000° T  +  0.002525°T2  +  0.000002T3  ,  (2.36b) 

p  =  334.328019°  +  4, 069. 032206° T  -  0.010344°T2  -  0.000012°T3  ,  (2.36c) 

where  T  is  the  number  of  Julian  centuries  (of  36,525  days)  reckoned  from 
January  0.5,  1900  at  Greenwich,  i.e.,  from  December  31,  1899,  12h  UT.  For 
January  0.0,  1978  at  Greenwich,  the  value  of  T  is  28,488.5/36,525;  upon 
considering  (2.36  a-c)  one  has 

[h]  =  279.310976°,  [s]  =  166.218322°,  [pj  =  268.055437°,  .  (2.37) 

where  the  brackets  have  been  used  to  indicate  this  specific  time  epoch.  Near 
a  point  of  expansion,  i.e.,  certainly  within  a  year,  h,  s  and  p  can  be  con¬ 
sidered  as  linear  functions  of  time  and  their  speeds  in  °/day,  etc.,  can  be 
evaluated  using  the  terms  linear  in  T  in  (2.36  a-c).  When  considered  together 
with  (2.37),  these  speeds  make  it  possible  to  compute  h,  s  and  p  for  any 
instant  in  1978  accurately  as 

h  =  279.310976°+  0.985647335°- D  +  0.04106864°-  hr 

+  0.00068448° «  min  +  0.00001141° -  sec  ,  (2.30a) 

s  =  166.218322°+  13. 176396769°- D  +  0.54901653° «  hr 

+  0.00915028°-  min  +  0.00015250°  -  sec  ,  (2.38b) 

p  =  268.055437°+  0. 111404030°.  D  +  0.00  464183°.  hr 

+  0.00007736°  .  min  +  0.00000129°  .  sec  ,  (2.38c) 

where  0  =  day  number  in  1978,  and  hr,  min,  sec  represent  hours,  minutes, 
seconds  in  UT  for  that  day.  From  the  formulas  (2.38  a-c)  the  various  rates 
of  change  in  h,  s  and  p  are  apparent.  They  also  confirm  the  periodicity  of 
h  (365.2421988  days  =  tropical  year),  of  s  (27.32158164  days  =  tropical  month) 
and  of  p  (8.847313  years  of  365.25  days). 


The  numerical  values  of  ctj  at  Greenwich  for  any  instant  in  1978 
can  be  found  from  the  general  expression  appearing  in  the  second  column  of 
Table  2.  The  rate  of  change  in  UT,  taken  in  the  interval  0-24  hours,  is 
15° /hr,  0.25° /min  and  0.00416667°/sec,  while  the  initial  values  and  the 
rates  of  change  in  the  other  three  variables,  h,  s  and  p,  have  been  given 
in  (2.38  a-c).  The  required  combinations  of  UT,  h,  s  and  p  thus  yield 


o^f  =  332.436645° 

ct'  =  258.162886° 
Mm 

aSSa=  iq8-621952° 


+  26.352793539°  *  D  + 
+  0.01830055°  «  min 

+  13.064992739°  *  D  + 
+  0.00907291°  *  min 

+  1.971294671°-  D  + 
+  0.00136895°  *  min 


1.09803306°- hr 
+  0.00030501°.  sec  , 

0.54437470°.  hr 
+  0.00015122°.  sec  , 

0.08213728°- hr 
+  0.00002282°- sec  ; 


(2.39a) 

(2.39b) 

(2.39c) 


cx^  =  9.310976°+  0.985647335°  .  D  +  15.04106864°.  hr 

1  +  0.25068448°  .  min  +  0.00417807°.  sec  , 

oig  =  216.874331°  -  25. 367146203°  «  D  +  13.94303557°- hr 
1  +  0.23238393°  -  min  +  0.00387307°- sec  , 

dp  =  350.689024°  -  0. 985647335°  »  D  +  14. 95893136°- hr 
1  +  0.24931552°  -  min  +  0.00415526°- sec  ; 


(2.40a) 

(2.40b) 

(2.40c) 


a'  =  226.185307°  -  24.  381498868°  »  D  +  28.98410421°. hr 
M? 

+  0.48306840°  «  min  +  0.00805114°«sec  ,  (2.41a) 

a'  =  30. °«  hr  +  0.5°«  min  +  0.00833333°  «  sec  ,  (2.41b) 

^2 
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=  328.022421°  -  37.446491607°  *  D  +  28.43972952°  «  hr 
2  +  0.47399549°  *  min  +  0.00789992°  «  sec  ,  (2.41c) 

a£  =  198.621952°  +  1.971294671°  *  D  +  30.08213728°  «  hr 

2  +  0.50136895°  «  min  +  0.00835615°  *  sec  ,  (2.41d) 

where  D,  hr,  min,  sec  were  defined  following  (2.38  a-c).  From  these  formulas 

the  rates  of  change  in  the  arguments  a',  and  thus  also  a.  are  apparent  and 

J  3 

agree  with  Table  2  of  [US]  wherever  they  are  comparable  (i.e.,  they  agree 
with  the  values  printed  in  [US]  as  "speed  per  solar  hour"  which,  however, 
exhibit  fewer  significant  digits  than  the  speeds  derived  above).  Further¬ 
more,  these  rates  also  agree  with  [Estes,  1980],  page  118,  and  with  [Godin, 
1972],  page  232;  they  agree  approximately  with  [Schwiderski ,  1980],  page  172, 
[Estes,  1980],  page  101,  and  [Lisitzin,  1974],  page  12.  The  rates  associated 
with  "D"  and  "hr"  are  further  presented  in  Table  2,  columns  3  and  4,  respective¬ 
ly,  under  the  headings  °/day  and  °/hour. 

The  fifth  column  of  Table  2  lists  a',  at  Greenwich  for  January  0.5, 

J 

1900  obtained,  with  the  aid  of  the  second  column,  from  (2.36  a-c)  for  T=0. 

One  could  evaluate  a',  at  any  instant  also  from  these  values  upon  applying  the 

J 

rates  listed  in  the  columns  3  and  4.  However,  this  would  lead  to  a  slight 
loss  of  accuracy  even  if  sufficient  digits  are  used  in  the  arithmetic,  due  to 
neglecting  the  terms  in  T2  (and  T3)  inherent  in  the  formulas  for  h,  s  and  p 
in  (2.36  a-c).  By  comparison,  the  terms  in  T2  and  T3  did  enter  (2.38  a-c) 
and  thus  also  the  expressions  that  followed,  developed  herein  in  view  of  SEASAT 
altimetry.  The  latter  formulas  are  advantageous  to  use  not  only  for  their 
accuracy,  but  also  because  they  are  very  simple  and  do  not  necessitate  a 
large  number  of  significant  digits  for  their  evaluation. 
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Having  found  all  of  the  aj  ,  one  can  recapitulate,  with  the  aid 
of  Table  2,  Greenwich  and  local  arguments  corresponding  to  SEASAT  observa¬ 
tional  epoch.  For  the  long-period  tides  there  is  no  distinction  between 


these  two  kinds 

and  we  have 

P 

o 

ii 

0  , 

(2.42a) 

aMf  = 

aMf  "  » 

(2.42b) 

aMm 

aMm  » 

(2.42c) 

aSSa  = 

aSSa  • 

(2.42d) 

Greenwich  arguments  for  the  diurnal  constituents  similarly  are 

a'  -  0.37°  , 

K1 

(2.43a) 

a'  +  0.50°  , 

U1 

(2.43b) 

ap  > 

(2.43c) 

while  Greenwich 

arguments  for  the  semidiurnal  constituents  are 

aM  = 

“m2  -  °-07°  • 

(2.44a) 

«S  = 
^2 

Otc  , 
b2 

(2.44b) 

aN  = 

<  -  0.07°  , 

^2 

(2.44c) 

aK  = 

c 

a'  -  0.66°  . 

K2 

(2.44d) 

The  above  three  groups  of  formulas  appeared,  in  a  similar  form,  in  Sections 
2.1,  2.2  and  2.3,  respectively. 
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The  local  arguments  (for  the  diurnal  and  semidiurnal  constituents) 


then  read 


ct’  +  X  -  0.37°  , 

h 

ciq  +  X  +  0.50°  , 


aA  +  X  ; 

K1 

+  2A  -  0.07°  , 


aA  +  2X  , 
z 

a'  +  2X  -  0.07°  , 
N2 

a'v  +  2X  -  0.66°  . 
*2 


(2.45a) 
(2.45b) 
(2.45c) 
(2.45d) 
(2.45e) 
(2.45f ) 
(2.45g) 
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3.  TIDAL  MATHEMATICAL  MODEL 


3. 1  Simplified  Tidal  Model 

This  model  encompasses  the  equilibrium  tide  in  conjunction  with 
the  solid  earth  deformation  as  described  on  pages  31  and  32  of  [Blaha,  198l] . 
In  that  report,  abbreviated  here  as  [B] ,  the  water  level  was  assumed  to 
conform  instantaneously  to  the  total  tidal  potential  f.  This  potential,  in 
turn,  was  assumed  to  be 

r  =  (l  +  k)r  ,  (3.1) 

where  r  is  the  equilibrium,  or  astronomical,  tidal  potential  (due  to  the 
effect  of  the  moon  and  the  sun)  and  k  is  one  of  the  Love  numbers;  kr  is  then 
the  "additional  tidal  potential"  owing  its  existence  to  the  earth  deformation. 
Equation  (3.1)  was  not  written  explicitly  per  se,  but  it  was  implied  in  the 
form  of  the  "geocentric  tide"  listed  below. 

The  simplifications  present  in  the  model  of  [B]  discarded  the 
forces  of  friction  and  viscosity  as  well  as  various  other  effects;  of  these, 
the  self-gravitation  and  the  ocean  loadinq  effects  will  be  discussed  in  the 
first  part  of  the  next  section.  The  heights  of  the  various  tidal  phenomena 
as  described  in  [B]  are  now  briefly  recapitulated: 

equilibrium  tide  =  17 g  ,  (3.2) 

earth  deformation  =  hr/g  ,  (3.3) 

additional  tide  (due  to  the  earth  deformation)  =  kr/g  ;  (3.4) 
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"geocentric  tide"  =  (l  +  k)r/g  ,  (3.5) 

"measured  tide"  =  (l  +  k-h)r/g  ;  (3.6) 

"geocentric  tide"  =  [(1  +  k)/(l  +  k  -  h)]  "measured  tide"  .  (3.7) 

The  symbol  g  in  these  formulas  represents  the  average  value  of  gravity. 
Equation  (3.7)  is  the  consequence  of  (3.5)  and  (3.6).  In  [B] ,  the  values 
of  the  Love  numbers  h,  k  were  adopted  from  [Vanfcek,  1980],  pages  10  and  11, 
as 

h  =  0.62  ,  k  =  0.29  .  (3.8) 
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3.2  Realistic  Tidal  Model 

In  this  model,  the  above  mentioned  forces  of  friction  and  viscosity 
as  well  as  the  self-gravitation  and  the  ocean  loading  effects  will  be  ac¬ 
counted  for  in  the  Laplace  Tidal  Equations  which  include  also  the  Coriolis 
force  and  the  spatial  distribution  of  depth.  The  development  will  begin  with 
the  self-gravitation  to  which  the  ocean  loading  effects  (the  effect  on  the 
ocean  bottom,  the  effect  on  the  potential)  will  be  related  subsequently.  The 
Laplace  Tidal  Equations  (LTE)  will  be  recapitulated  using  a  fair  amount  of 
detail,  which  will  shed  light  on  various  possible  simplifications  and,  especial¬ 
ly,  on  the  surface  tide  directly  related  to  altimeter  observations.  Finally, 
a  model  expressing  the  surface  tide  for  the  diurnal  and  semidiurnal  tidal 
constituents  of  interest  will  be  derived.  The  following  four  references  will 
be  mentioned  on  numerous  occasions:  [Vam'cek,  1980],  [Estes,  1980],  [Parke 
and  Hendershott,  1980]  and  [Schwiderski ,  1980];  they  will  be  abbreviated  in  the 
text  as  [v] ,  [E],  [PH]  and  [S] ,  respectively. 

Effect  of  the  self-gravitation  of  tidal  waters  on  the  potential.  We 
begin  by  presenting  the  effect  of  (discontinuous)  point  masses  on  the  poten¬ 
tial  at  point  i,  denoted  Vs: 

Vs  =  G  I  . 

where 

G  =  gravitational  constant, 

M.  =  mass  of  the  j-th  point  mass, 

J 

2--  =  distance  between  the  point  i  (of  evaluation) 

J  and  the  j-th  point  mass. 
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If  is  replaced  by  PjCjdftj,  where  p  is  the  average  density  of  mass  (here 
water)  in  the  column  of  height  £  and  cross-section  dft,  the  potential  becomes 

ws  ■  G  \  <®j  .  (3.9a) 

Since  water  has  continuous  distribution  we  can  write 

V  =  G  JX  PC  |x-  x'l'1  dS'  ,  (3.9b) 

s 

where,  in  addition  to  the  symbols  already  introduced  (e.g.,  £  represents  the 
height  of  the  tidal  waters),  the  following  representations  apply: 

x  ...  position  where  V$  is  being  evaluated  , 
x1  ...  position  of  the  element  of  integration  , 
dS'...  surface  element  , 

S  ...  integration  domain,  the  ocean  surface. 

Equation  (3.9b)  is  essentially  (55)  of  [V]  (the  signs  in  this  reference  indicate 
"correction"  rather  than  "effect"). 

The  simplifications  to  be  used  throughout  are  the  spherical  approxima¬ 
tion  and  the  assumption  of  constant  density  of  the  sea  water.  The  following 
relations  thus  hold  true  only  approximately,  but  the  sign  of  equality  can 
nevertheless  be  used  if  one  keeps  in  mind  this  limitation.  Under  these  circum¬ 
stances,  from  (3.9b)  we  obtain 

V  •  GR2p  JJ  £  |  x  -  x 1 1  ” 1  dSV  ,  (3.10) 

fl 
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where 


R  =  mean  radius  of  the  earth  , 

p  =  mean  density  of  the  sea  water  , 

dfi1  =  solid  angle  element  (e.g.,  cos^'d^'dA' )  , 

=  ocean  surface  over  the  unit  sphere  . 

Equation  (3.10)  is,  in  fact,  the  expression  for  on  page  105  of  [E]  .  All  the 
relations  written  in  terms  of  G  and  R  could  also  be  written  in  terms  of  g  accord¬ 
ing  to 

g  =  GMe/R2,  (3.11a) 

or 

GpR  =  gR3p/Me  ,  (3.11b) 

where  fl  is  the  earth's  mass, 
e 

The  distance  7,  7'  can  be  expressed  as 

| x  -  x* |  =  R  j x  -  x' |  ,  (3.12) 

where  the  points  x,  x'  are  now  located  on  a  unit  sphere;  either  pair  is 
separated  by  the  central  angle  ^  such  that 

cos4>  =  sin4>  sin4>'  +  cos4>  cos<J>'  cos(A' -  A)  ,  (3.13) 

where  <p, A  are  the  latitude  and  the  longitude,  respectively.  The  reciprocal 
value  of  the  distance  x,  x1  may  be  expanded  in  an  infinite  series  of  Legendre 
polynomials  in  the  argument  cos^,  Pp(cos^),  so  that  we  can  write 

( 3. 14) 

|x- x '  j  ~ 1  =  (2  sin  4^)_1  =  [2(1- cos^)]'*5  =  l  P  (cos^)  . 

n=0 
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With  (3.12)  and  (3.14),  equation  (3.10)  becomes 

V$  »  (gR3p/Me)JJ  c(2  sin  Js*)"1  da'  =  (gR3p/Me) £  [2(1  -  cos.j,)]'1'2  d ft' 

ft  ft 

*  (gR3p/Me)/JC  I  Pn(cos^)  diT  ,  (3.15) 

ft  n 

where  (3.11b)  has  also  been  taken  into  account  and  where  the  notation 

co 

l  *  I 

n  n=0 

has  been  introduced.  The  first  expression  in  (3.15)  corresponds  to  a  formula 
on  page  106  and  to  one  building  block  of  (A10)  in  [E] ,  the  second  expression 
corresponds  to  (59)  of  [V] ,  and  the  third  expression  corresponds  to  (56)  of 

[V]- 

Because  the  integrand  in  the  last  expression  of  (3.15)  is  a  well- 
behaved  function,  the  summation  can  be  interchanged  with  the  integration  and 
we  have 

Vs  =  (gR*p/Me)  [ffe  Pn(cos^)  dir  =  l  (V  )  ,  (3.16a) 

n  ft  n 

(V$)n  =  (gR3o/Me)jy^p  (Cosi|>)  diT  .  (3.16b) 

ft  n 

In  these  equations  the  potential  V$  is  expanded  in  terms  of  surface  spherical 
harmonics  (this  is  always  possible  for  a  function  on  a  sphere)  and  its  explicit 
form  is  given. 

The  vertical  displacement  of  an  equipotential  surface  caused  by  the 
self-gravitation  is 

ua  =  V$/g  ,  (3.17) 
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which  is  the  same  as  equation  (60)  in  [V] .  Upon  dividing  (3.15)  by  g,  one 
obtains 

u.  «  (R3p/Mjjn  (2  sin  Hip)'1  dO'  ,  (3. 17' ) 

d  e  a 

which  could  be  further  expanded  to  include  the  other  two  expressions  in 
(3.15).  From  (3.16a,b)  it  similarly  follows  that 

ua  =  (R3P/Me)  l  tit  Pn(cos^)  dQ1  =  l  (ua)n  ,  (3.18a) 

n  ft  n 

(ua}n  =  (VsV9  =  (RVMe)£fS  Pn  (cos^)  dn‘  .  (3.18b) 

The  formula  giving  (u  )  can  be  rewritten  if  one  takes  advantaae  of 
the  result  from  Appendix  2  (see  A2.6)  derived  in  the  spherical  coordinates 
(9,A),  where  0  is  the  colatitude.  If  p£  further  denotes  the  average  earth 
density, 

Pe  =  Me/[(4/3)7TR3]  , 
equation  (3.18b)  becomes 

(ua)n  =  [3P/(47TPe^]-tf  C(O'.X')  Pn  (cosip)  dfi'  . 

But  according  to  Appendix  2  the  last  integral  is  [4tt/ ( 2n  +  1)] Cn ( © » A ) ,  hence 
(ua}n  =  t3/(2n+1)l  (P/Pe)  » 

where  £n  =  £n(e,A)  is  the  surface  spherical  harmonic  of  n-th  degree,  associated 
with  £.  We  can  thus  write 
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(3.19a) 


and 

“a  -  l  %  ?„  <3-19b> 

n 

where 

(3.19c) 

=  [3/(2n+  1)]  *  (average  water  density/average  earth  density). 

The  expression  for  an  agrees  with  [PH],  page  393,  if  one  disregards  an  apparent 
typographical  error  in  that  the  first  factor  therein  is  written  as  (3/2n+l). 

Effect  of  the  ocean  tidal  loading  on  the  ocean  bottom.  The  deforma¬ 
tion  of  the  ocean  bottom  due  to  the  ocean  tidal  loading  is  symbolized  by  u^. 

In  our  context,  it  is  a  function  on  a  sphere  which  can  be  expanded  in  surface 
spherical -harmonics  as 

"  l  <"A  '  <3-20> 

n 

As  is  explained  in  [v] ,  pages  21  and  22,  the  displacement  ua  serves  as  a  norm 
for  u^,  giving  rise  to  the  following  load  numbers  (sometimes  called  Love  load 
numbers)  for  all  possible  wave  numbers  n: 


This  together  with  (3.18b)  and  (3.20)  yield 

(3.21) 

utsI  hn(ua)n  =  l  hn(Vn  =  ^R3p/Me^  E  hn#£pn  (cos^  dfi'  » 

n  n  n  Q 


where  the  second  equality  corresponds  to  (64)  in  [V].  If,  in  the  last  expres¬ 
sion,  the  summation  and  the  integration  are  interchanged  one  has 
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(3.22a) 


u£  =  (R3P/Me)/J"5l  hn  pn  ( cos^j)  dfi*  =  R2pXf  £  U1  (i|0  dfi'  , 

fin  fi 


U'(i|»)  =  (R/Me)  l  h‘  Pn  (cosip)  . 

n 


(3.22b) 


The  first  equality  in  (3.22a)  corresponds  to  (65)  and  (66)  of  [v] *  the  second 
equality  serves  as  another  building  block  in  (A10)  of  [E] ,  and  (3.22b)  is  the 
definition  appearing  in  (A5)  of  [E] . 

If  (u  )  in  the  first  equivalence  of  (3.21)  is  taken  from  (3.19a), 

a  n 

we  write  at  once: 


ul  ^  hn  an  ?n  ’ 


(3.23) 


which  corresponds  to  the  “ocean  induced  vertical  component  of  the  solid  earth 
tide"  as  it  appears  on  page  393  of  [PH]. 

Effect  of  the  ocean  tidal  loading  on  the  potential.  The  disturbing 
potential  due  to  the  deformation  u^  just  treated  is  denoted  as  Vb.  Similar 
to  (3.17),  it  gives  rise  to  the  displacement  u^  of  an  equipotential  surface, 


ui  ■  v«  ■ 


(3.24) 


This  displacement  can  again  be  expanded  in  surface  spherical -harmonics: 


ui  =  l  <“i>n  • 


(3.25) 


The  statement  that  followed  (3,20)  now  applies  to  u.  in  conjunction  with 
another  kind  of  load  numbers: 


K  -  (ui)r./(ua>n  • 
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This  together  with  (3.18b)  and  (3.25)  yield 


(3.26) 


«i  =  l  kn  (ua)n  =  (l/g>  l  k^  (Vs)n  =  (R3p/Me)  l  Pn  (cos*)  dft1  , 


where  the  second  equality  corresponds  to  (64)  of  [V] .  Upon  interchanging  the 
summation  and  the  integration  in  the  last  expression  we  obtain 


u.  =  (R3p/Me)  Jjf  5  l  k;  Pn  (cos*)  dn*  =  (R2p/g)  XT  K  *'(*)  dfl'  , 

fin  fi 

<&'(*)  =  (Rg/Me)  I  k^  Pn  (cos*)  . 


(3.27a) 

(3.27b) 


The  last  expression  in  (3.27a)  serves  as  a  third  building  block  in  (A10)  of 
[E] ,  and  (3.27b)  is  the  same  definition  as  in  (A5)  of  [E] .  From  (3.24)  and 
(3.27a)  we  also  have 

Vb  =  gu.  =  R2p/A  (4i)  da1  •  (3.28) 

Q 

If  (u  )  in  the  first  equivalence  of  (3.26)  is  taken  from  (3.19a), 
'an 

it  follows  that 


u 


i 


(3.29) 


Upon  utilizing  this  form  for  ,  one  obtains 

\  ■  3  l  K  “n  5„  .  (3-3°) 

n 

which  corresponds  to  the  second  part  of  the  "ocean  induced  potential  at  the 
mean  sea  surface"  as  presented  on  page  393  of  [PH] . 
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The  total  tidal  potential  causing  the  water  to  depart  from  its 
average  level  is  composed  of  the  astronomical  tidal  potential  r,  of  the 
additional  tidal  potential  kr  due  to  the  earth  deformation,  of  V$  due  to 
the  self-gravitation,  and  of  due  to  the  ocean  loading: 

‘tot  ■  <1+k>r  ♦  vs  +  vb  •  <3-31> 

When  written  in  terms  of  (3.17)  and  (3.24),  this  becomes 

rtot  =  (1+  k)r  +  gua  +  gui  ,  (3.31') 

which,  with  (3.19b)  and  (3.29),  also  is 

rtot  ■  9  1  ,ltkn>  “n  • 

n 

this  formulation  corresponds  to  that  on  page  393  of  [PH] .  If  the  water  level 
adjusted  itself  instantaneously  to  an  equipotential  surface  implied  by  (3.31'), 
its  height  above  the  mean  surface  would  be 

(3.32) 

This  quantity  would  not  be  the  same  as  £  obtained  from  the  solution  of  the 
ITE,  involving  the  forces  of  friction,  viscosity,  etc.  This  fact  is  reflected 
by  the  prime  attributed  to  in  (3.32). 

The  departure  of  the  ocean  bottom  from  its  average  value  is  given  as 
£b  =  hr/ g  +  u£  .  (3.33) 
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If  u^  is  written  as  in  (3.23),  equation  (3.33)  becomes 
«b  -  hr/9  +  l  h;  “n  ■ 

n 

which  again  corresponds  to  the  formulation  on  page  393  of  [PH] .  Under  the 
simplification  (3.32),  the  ocean  tide  would  be 

V  =  V$  -  q  *  (l+k-h)r/g  +  (ua+u.  -Ujl).  (3.34) 

If,  in  addition,  we  also  ignored  the  self-gravitation  and  the  ocean 
loading  effects,  i.e.,  u  ,  u.  and  u.  in  (3.32)-(3.34),we  would  have 

^  3  1  jL 

q  =  (1 +  k)r/g  , 

q  =  hr/g  , 

r  =  q  -  q  =  d+k-h)r/g  . 

But  this  is  the  simplified  tidal  model  used  in  [B]  and  recapitulated  in  the 
previous  section,  where  the  earth  deformation  corresponds  to  the  above  ££  , 
the  "geocentric  tide"  corresponds  to  q  ,  and  the  "measured  tide"  corresponds 
to  £"  .  The  above  process  of  successively  greater  simplifications  has  thus 
brought  into  focus  the  assumptions  and  simplifications  present  in  the  model 
of  [B]. 

Most  of  the  quantities  just  discussed  are  shown  in  Figure  1.  The 
figure  represents  the  realistic  tidal  model  and  involves  the  solution  of  the 
LTE;  this  fact  is  illustrated  by  the  presence  of  the  quantity  y.  A  brief 
description  of  the  quantities  in  Figure  1  follows: 
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r/g  =  equilibrium  tide  , 

y  =  quantity  resulting  from  the  solution  of  the  LTE  , 

kr/g  =  additional  tide  (due  to  the  earth  deformation  hr/g)  , 

u-  =  V./g  =  displacement  due  to  the  ocean  tidal  loading  on 

1  the  ocean  bottom  (i.e.,  due  to  u^)  , 

u  =  V  / g  =  displacement  due  to  the  self-gravitation  of 

a  s  tidal  waters  ; 

hr/g  =  earth  deformation  (called  also  solid  earth  tide)  , 
u^  =  ocean  bottom  deformation  due  to  the  ocean  tidal  loading  ; 

£  =  surface  tide  , 

s 

£b  =  bottom  tide  , 

C  =  ocean  tide  . 


From  Figure  1  it  follows  that 

£s  =  (1  +  k)r/g  +  y  +  ua  +  u.  ,  (3.35) 

=  hr/g  +  u^»  (3.36) 

£  *  Cs  -  Cb  =  (1  +  k-  h)r/g  +  y  +  (ug  +  ui  -  u4)  .  (3.37) 

It  is  also  apparent  from  the  figure  that 

£  =  h(instant.)  -  h(ave.)  .  (3. 37') 


Without  the  presence  of  y,  (3.35)-(3. 37)  would  reduce  to  (3.32-  3.34). 
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ocean  tide) 


At  this  point,  three  remarks  can  be  made  upon  examining  Figure  1. 

First,  identifying  the  average  water  surface  with  the  geoid  indicates  a  few 
approximations  (such  as  considering  the  average  water  surface  to  be  equi- 
potential,  which  is  not  exactly  fulfilled  in  several  oceanic  areas)  and,  more 
importantly,  it  indicates  that  the  permanent  tide  is  ignored  in  the  present 
considerations.  However,  this  tidal  effect,  which  essentially  transforms 
the  "static  geoid"  (of  interest)  into  the  above  geoid,  is  treated  separately. 

In  fact,  all  the  pertinent  long-period  tidal  effects  are  envisioned  to  be 
treated  separately  along  the  lines  of  the  simplified  model;  the  realistic 
tidal  model  as  adapted  from  [E]  was  developed  only  for  the  principal  diurnal 
and  semidiurnal  tides.  We  note  that  in  [S]  ,  pages  165  and  166,  and  in  [Estes, 
1977],  Figure  1,  the  same  identification  (ave.  water  surface  =  geoid)  was 
made  as  in  our  Figure  1. 

Second,  we  remark  that  for  £>0,  one  must  have  u^<0  and  thus  also 
u^<0;  the  response  of  u^  to  u£  is  qualitatively  somewhat  similar  to  the  re¬ 
sponse  of  the  additional  tide  (kr/g)  to  the  earth  deformation  (hr/g).  On  the 
other  hand,  if  £>0,  the  disturbing  potential  increases  (see  e.g.  the  illustra¬ 
tive  case  3.9a)  and  thus  the  corresponding  geoid  undulation  also  increases, 
i.e.,  one  has  u  >0.  Therefore,  u.  and  u„  have  the  opposite  signs  from  that  of 
u,  and  should  have  been  directed  downwards  in  Figure  1  where  £  and  u,  have  an 

a  a 

upward  sense.  In  Figure  1  of  [S] ,  u^  is  indeed  represented  in  this  manner  and 
is  appropriately  called  the  solid  earth  dip,  £eo  (i.e.,  the  earth-dip  response 
to  the  oceanic  tidal  load  £).  But  this  is  merely  a  matter  of  convention  since 
in  this  reference  £b  is  expressed  essentially  as  hr/g  -  (positive  quantity,  £e0) 
while  in  the  present  study  is  given  as  hr/g  +  (negative  quantity,  u^).  A 
similar  comment  would  also  apply  to  u. ,  as  may  be  gathered  upon  consulting 
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equation  (9)  of  the  same  reference.  Figure  10  of  [Vl  where  k'  and  h'  are 
seen  to  be  negative  for  all  values  of  n  helps  us  to  verify,  in  a  plausible 
way,  that  the  quantities  ui  and  u^  have  the  opposite  signs  of  ug  . 

The  last  remark  concerns  the  notations  used  in  the  figure  as  well 
as  throughout  the  text.  As  stated  in  [V],  page  10,  only  the  surface  values 
h2  ,  k 2  of  Love  numbers  hn  ,  kn  are  used  in  geodetic  applications,  in  which 
case  the  symbols  h  ,  k  are  often  employed.  These  symbols  have  been  adopted 
in  the  present  study.  It  was  demonstrated  on  pages  3  and  4  of  [V]  that  the 
astronomical  tidal  potential,  r,  can  for  most  practical  purposes  be  equated 
with  r2  of  the  development  in  terms  of  Legendre  polynomials.  It  should  thus 
be  borne  in  mind  that  in  our  context  r,  kr  and  hr  could  replace  r2  ,  k^ 
and  h2r2  or  l  r  ,  £  kp  Tn  and  £  hn  Tn  ,  respectively,  used  in  some  references. 

n  n  n 

In  the  upcoming  derivation,  the  horizontal  components  of  the  force 
corresponding  to  (3.31)  or  (3.31')  will  be  needed  along  the  X-  and  <j>- coordinate 
lines.  Since  one  can  write  for  this  force: 

T'  =  9rad  rtot  . 

the  required  components  are 

(3.38a) 

F\=  [1/(R  C0S4>)]3rt0t/3A  «  [1/(R  cosc}))][(l+k)9r/3A+g3(ua+ui)/9X]  , 

(3.38b) 

=  ( l/R)  artot/ad>  =  (l/R)  [(l  +  k)  ar/a<j>  +  ga(ua  +  ui )/B4>] . 

These  components  will  be  used ,  together  wi th  the  components  F^"  and  F^'  of  F ", 
the  force  of  friction  and  viscosity,  in  the  LTE  in  the  form  (Al)  of  [E] .  The 
components  of  F"  will  be  transcribed  from  this  reference,  but  otherwise  the 


final  form  of  these  equations  will  be  reden' ved  independently,  for  the 
purpose  of  scrutinizing  the  quantity  £s  needed  for  a  refinement  of  a  1 ti - 
metric  applications.  The  model  of  [E]  is  advantageous  especially  because 
the  tidal  elevations  £  for  all  the  diurnal  and  semidiurnal  constituents  of 
interest  are  developed  in  terms  of  spherical -harmonic  tidal  coefficients, 
tailored  for  use  in  the  first  phase  of  satellite  altimetry  adjustment 
(global  adjustment)  with  the  spherical-harmonic  potential  coefficients  as 
parameters . 

The  first  two  (out  of  three)  LTE  presented  in  (Al)  of  [E]  read 
3u/3t  -  2oj  sin^i  v  =  -[g/(R  cos<}>)] 35$/3A  +  +  F^'  ,  (3.39a) 

3v/3t  +  2u)  si n<p  u  =  -  (g/R)3£s/3<t>  +  F^  +  F^'  ,  (3.39b) 

where,  in  addition  to  the  symbols  already  defined,  u  and  v  are  the  eastward 
and  northward  components  of  the  water  velocity  and  to  is  the  rotational  veloc¬ 
ity  of  the  earth.  These  equations  involve  not  only  the  tidal  generating  po¬ 
tential  (it  gives  rise  to  ,  F^') ,  but  also  the  forces  of  friction  and 
viscosity  (see  F^'  ,  F'^  mentioned  earlier)  and  the  Coriolis  force  (it  is 
associated  with  to  on  the  left-hand  sides  of  3.39  a,b).  The  third  equation 
in  (Al)  involves  the  water  depth  and  it  will  be  transcribed  directly  into  the 
final  formulas.  The  three  equations  (Al)  correspond  to  equations  (5)  -  (7) 
of  [Lisitzin,  1974]. 

Next,  F^'  and  F  1  from  (3.38  a,b)  are  substituted  into  (3.39  a,b), 
together  with  £$,  expressed  from  Figure  1  as 

=  K  +  5b  *  £  +  hr/g  +  \x%  .  (3.40) 
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This  results  in 


(3.41a) 


3u/3t  -  2w  sin*  v  =  [l/(R  cos*)]  [-g3C/3A  +  (1  +  k-h)  3T/3A 

+  q8(ua  +  u.  -  u£)/3A]  +  Fx”  , 

(3.41b) 

3v/3t  +  2ui  sin*  u  =  (l/R)[-g3?/3*  +  (l  +  k-h)3r/3* 

+  g3(ua  +  -  u^)/3*]  +  FJ'  . 

According  to  (3.19b),  (3.23)  and  (3.29),  we  could  write 
ua  +  ui  '  ujl  =  I  (ltk;-|’n)“r,5n  • 

n 

However,  to  relate  this  quantity  to  the  results  in  [E] ,  we  use  (3. 17'),  (3.22a) 
and  (3.27a): 

-i  ,  (3.42) 

ua  +  ui  ■  =  R2p  $  ?  [(R/M  ) (2  sin  *#)  +  (1/g)*'  (*)  -  U'  (*)]dft‘ . 

3  1  £1  e 

In  agreement  with  [E] ,  I”  is  defined  at  (*,A,t)  as 

T1  =  g(ufl  +  u.  -  u£)  .  (3.43a) 

With  the  aid  of  (3.42),  equation  (3.43a)  reads 

T'  »  .f/C"  "Green"  d£I'  ,  (3.43b) 

Q 

where  the  Green's  function  is 

"Green"  =  R2p[Rg/Mg)  (2  sin  ^*)_1 +  $'(*)- qU' (*)]  .  (3.43c) 

The  last  two  equations  correspond  to  the  result  on  pages  106  and  107  of  [E]  , 
especially  equation  (A10). 


The  last  step  consists  of  utilizing  (3.43a)  in  (3.41  a,b)  and 


transcribing  F^1  ,  F^1  from  [E] ,  page  107;  the  two  equations  are  then  comple¬ 
mented  by  the  third  equation  in  (Al),  page  93  of  [E] .  The  result  is 

9u/3t  -  2w  sin<}>  v  =  [1/(R  cos<i>}]  [-q3C/3X  +  (1  +  k  -  h)  3r/3X+9r'/9A] 

-  Cr(u2  +  v2)}$  u/H  +  ChvAu  ,  (3.44a) 

3v/3t  +  2w  sin4>  u  =  ( 1/R)  [  -  g3£/3<j>  +  ( 1  +  k  -  h)  3r/S<J>  +  3T ' /3d>] 

-  Cr(u2  +  v2)'2  v/H  +  ChvAv  ,  (3.44b) 

3£/3t  =  - [l/( R  cos <J>) J  [3(Hu)/3X  + 3(H  vcos4>)/3(j)]  ,  (3.44c) 

where  A  is  the  horizontal  Laplacian,  Cr  and  Chv  are  the  coefficients  of  friction 
and  horizontal  eddy  viscosity,  and 

H  =  h(ave.)  +  £  =  h(instant.)  . 

Detailed  information  about  the  numerical  solution  of  equations  (3.44  a-c)  can 
be  found  on  pages  106-117  of  [E] . 

After  C  has  been  solved  for,  the  quantity  directly  related  to 
satellite  altimetry,  can  be  expressed  as  in  (3.40),  namely 

=  5  +  hr/g  +  u£  ,  (3.45a) 

where 

hr/g  =  h  »  (equilibrium  tide)  .  (3.45b) 

In  practical  adjustments  of  altimeter  data,  the  quantity  u^  is  envisioned 
either  to  be  neglected  or  to  be  expressed  very  approximately  as  some  multiple 


of  £.  The  latter  choice  could  be  carried  out  in  conjunction  with  the  M2 
constituent  for  which  £  in  [E]  is  expressed  through  the  solution  of  (3.44  a-c) 
in  an  iterative  process  (in  this  case  r  is  made  to  correspond  to  M2).  The 
remaining  semidiurnal  and  diurnal  constituents  are  treated,  in  this  refer¬ 
ence,  in  a  more  approximate  fashion  by  ignoring  r'.  In  order  to  express  £$ 
as  accurately  as  practicable  the  quantities  u,  ,  u.  and  u0  as  well  as  their 
combinations  (u  +u. )  and,  especially,  u  +u. -u  =r'/g  from  (3.43a)  will  be 

a  I  a  1  x, 

examined  in  the  light  of  the  numerical  solutions  presented  in  [PH]  and  [E] . 

The  contour  maps  in  Figures  1,  3  and  5  of  [PH]  depict  £  for  the  con- 

and  K,,  respectively.  The  contour  maps  of  (u  +u.)  for  these 

1  a  1 

constituents  are  shown  in  Figures  10,  11  and  12,  and  the  contour  maps  of  u^ 
are  similarly  shown  in  Figures  13,  14  and  15  of  this  reference.  When  examin¬ 
ing  the  results  for  we  imagine  the  phase  changed  by  180°  and,  accordingly, 
the  sign  changed.  Upon  inspection,  the  quantity  -u^  is  then  seen  to  be  es¬ 
sentially  "in  phase"  with  £,  consistent  with  our  earlier  discussion. 

For  the  M2  constituent,  €  ranges  mostly  within  25  -  50  cm  (exceptionally, 
it  reaches  over  lm),  while  (u  +u. )  ranges  mostly  within  1-3  cm  (exceptionally 

a  1 

4  cm),  and  -u^  ranges  mostly  within  1-4  cm  (in  two  coastal  areas  it  reaches 

about  6cm).  The  phases  as  well  as  the  highs  and  lows  of  these  quantities  have 

similar  characteristics  over  most  of  the  global  oceans.  From  this  inspection 
one  can  draw  very  approximate  tentative  conclusions: 

u„  +  u.  <  0.05£,  -u„  >0.05£;  u,  +  u.  -  u„  =  0.10£  .  (3.46) 

Nearly  the  same  conclusions  can  be  reached  also  fo>"  the  S2  and  Kj  constituents. 
In  the  case  of  S2  ,  £  is  mostly  in  the  10  -  30  cm  range  (exceptionally  40  cm), 


stituents  M2,  S2 
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( Ua  +  ui )  1S  in  t^ie  -1.5  cm  range  and  -u^  is  mostly  in  the  0.5 -1.5  cm 
range  (exceptionally  2cm).  The  ranges  for  the  constituent  are  similar  to 
these. 

Next,  [E]  is  consulted  with  regard  to  the  constituent  for  which 
contour  maps  are  presented.  Figures  A1  (for  £),  A8  (for  u  +u--u0),  A9  (for 

d  1  36 

u  )  and  A10  (for  -u„)  allow  the  following  observations:  £  is  mostly  in  the 
10  -  50  cm  range  (exceptionally  80  cm  in  the  Equatorial  Pacific  and  east  of 
Australia);  u  +  u.  -  u„  is  mostly  in  the  1  -  5  cm  range  (exceptionally  12  cm  in 
the  Equatorial  Pacific);  and  u,  and  u„,  both  having  very  similar  characteris- 
tics,  are  typically  in  the  1  -  3  cm  range  (exceptionally  7  cm  in  the  Equatorial 
Pacific).  There  is  a  generally  good  agreement  in  the  phases  and  in  the  rela¬ 
tive  magnitudes  of  these  four  quantities.  A  visual  inspection  leads  to  the 
following  approximate  conclusions: 

u  > 0.05^  ,  -u„  >0.055  ;  u,  +  u.  -  J0-  0.105  .  (3.47) 

a  30  a  I  36 

The  last  relations  in  (3.46)  and  (3.47)  are  supported  by  the  result  described 
on  pages  170  and  171  of  [S].  Transcribed  into  our  notations  it  states  that 

ua  +  u.  -  u£  =  r'/g  «  0.105  , 

obtained  "after  evaluating  the  Green's  function  representation  of  the  three 
oceanic  tidal  load  effects"  (page  171).  On  the  same  page  of  [S]  it  is  also 
stated:  "The  author's  computations  supported  the  marginal  effects  of  oceanic 
tidal  loading  found  by  Estes  (1977)." 

From  the  above,  especially  from  the  contour  maps  associated  with  the 
constituent,  the  various  effects  can  be  estimated  very  approximately  as 


follows: 


u  -  0.06C  ,  u.  =  -0.02C  ,  u  +u.  »  0.04C  ; 

a  I  3  1 

-uA»  0.06C  ;  (3.48) 

u.  +  u.  -  u0  =  0. 10C  . 

a  I 

Of  these,  (3.48)  is  the  most  important,  allowing  us  to  express  C$  from 
(3.45  a,b)  as 

Cs  =  0.94C  +  h  *  (equilibrium  tide)  .  (3.49) 

This  formula  could  be  used  in  conjunction  with  any  constituent  whose  C  is  com¬ 
puted  rigorously  as  in  (3.44  a-c).  In  [e]  such  a  computation  was  performed 
for  the  constituent.  As  was  already  mentioned,  the  self-gravitation  and 
the  ocean  loading  effects  were  neglected  in  [E]  (I”  was  disregarded)  during 
the  computations  of  C  for  the  remaining  diurnal  and  semidiurnal  constituents. 

To  be  consistent  with  this  procedure  when  adopting  the  spherical-harmonic  tidal 
coefficients  from  [E] ,  we  neglect  u^  in  (3.45a)  applied  to  these  constituents, 
in  which  case  (3.45  a,b)  yield 

C$  =  T  +  h  *  (equilibrium  tide)  .  (3.50) 

The  above  analysis  suggests  the  use  of  (3.49)  in  conjunction  with 
M2  and  the  use  of  (3.50)  in  conjunction  with  all  othe  other  diurnal  and  semi¬ 
diurnal  constituents  considered.  In  fact,  (3.50)  could  have  been  designed 
for  M2  as  well  since,  according  to  the  earlier  results  taken  from  [PH]  and 
[E] ,  the  magnitude  of  u^  is  usually  well  below  the  5cm  level  and  the  effect 
of  u^  could  also  qualify  as  being  only  marginal.  However,  we  prefer  to  use 
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a  more  accurate  model  (3.49)  because  not  only  the  tidal  amplitudes,  but  also 
the  phase  angles  are  now  considered  to  be  subject  to  adjustment,  and  the  model 
is  nonlinear  in  the  phase  angles.  Hence  the  best  possible  initial  values 
should  be  sought  for  the  parameters  appearing  in  the  linearized  observation 
equations.  Although  one  could  argue  that  the  value  5  is  not  significantly 
different  from  the  value  0.945  for  this  purpose,  it  is  equally  true  that  the 
more  accurate  model  for  represented  by  (3.49)  is  as  easy  to  handle  as  the 
model  represented  by  (3.50). 

As  was  done  by  Estes  [1980]  for  the  desired  diurnal  and  semidiurnal 
constituents,  the  ocean  tide  5  can  be  expressed  by  means  of  spherical -harmonic 
tidal  coefficients  encompassing,  through  a  least-squares  fit,  the  solution  of 
the  LTE.  For  a  given  tidal  constituent,  these  coefficients  lead  to  an  expansion 
similar  to  equation  (3)  of  [Estes,  1980].  This  formulation  is  described  with 
the  aid  of  the  following  model: 

5.  =  A.  cos(a. -tjy.)  =  A.  cos*];,  cosct.  +  A.  sinO.  sina.  ,  (3.51) 

J  J  JJ  J  J  J  J  vJ  J 

where 

5j  =  constituent  height  observable  by  tidal  gauges, 
a.  =  Greenwich  argument  of  the  constituent, 

J 

A.  =  A.(<j>,A)  =  amplitude  of  the  constituent, 

3  3 

=  xp - ( <j>, X)  =  phase  angle  of  the  constituent. 

The  longitude  (A)  of  the  place  of  observation  is  not  explicitly  needed  since 

it  is  included  in  <p..  The  angles  a.  are  thus  computed  as  in  (2.43a)  -  (2. 44d). 

0  0 
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Equation  (3.51)  is  reformulated  to  read 


£•  =  a.  cosa.  +  b.  sina.  , 

J  J  J  J  J 

where 

a.  =  a.(<j>,X)  =  A.  cosij).  =  l  l  (a.  cos  mX  +  b.  sin  mX)P  (sin<}>), 

J  J  J  J  n  m  Jnm  Jnm 

b.  =  b .  (<J>,X)  =  A.  sinijj.  =  l  £  (c.  cos  mX  +  d.  sin  mX)P  (si n<f») , 

J  J  J  J  n  m  Jnm  Jnm  nm 


(3.52a) 

(3.52b) 

(3.52c) 


from  which  it  follows  that 

Aj  =  (a?  +  b?)*5  ,  (3.52d) 

cos^.  =  a  ../A,.  ,  sin^  =  b^/Aj  .  (3.52e) 

In  these  formulas  a^  ,  etc.,  are  the  spherical -harmonic  tidal  coefficients 

Jnm 

of  degree  and  order  (n,m)  belonging  to  the  constituent  j,  and  pnm(sin<l>)  are 
the  associated  Legendre  functions  in  the  argument  sin<f>,<)>  being  the  geocentric 
latitude.  In  the  next  chapter,  the  above  formulation  will  be  used  in  conjunc¬ 
tion  with  the  model  for  seen  previously. 
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4.  TIDAL  ADJUSTMENT 


4. 1  Inclusion  of  Tidal  Effects  in  the  Altimetry  Adjustment  Model 

The  basic  model  equation  of  satellite  altimetry  was  written  in 
equation  (3.1)  of  [Blaha,  1979]  as 

H  =  R  -  r  +  d  , 

where  H  represents  the  altimetry,  R  is  the  distance  from  the  geocenter  to 
the  satellite  at  the  time  of  observation,  d  is  a  correction,  always  smaller 
than  5m  for  the  satellite  altitude  under  1,000  km  as  described  e.q.  on 
page  28  of  [Blaha,  1977]  or  in  [Blaha,  1977'],  and  r  is  the  distance  from 
the  geocenter  to  a  sub-satellite  point  on  the  sea  surface;  it  is  given  on 
page  15  of  [Blaha,  1979]  as 

r  =  r'  +  N  , 

where  r1  is  the  corresponding  distance  to  the  (geocentric)  reference  ellipsoid 
and  N  represents  the  geoid  undulation.  The  main  feature  of  an  earlier  approach 
consisted  in  expressing  N  (and  thus  r)  in  terms  of  the  geoidal  parameters  only, 
as  if  the  measured  sea  surface  coincided  with  the  geoid.  Although  this  model 
deficiency  was  of  little  consequence  in  past  adjustments  of  GE0S-3  altimetry, 
it  will  be  at  least  partly  removed  from  the  SEASAT  altimetry  model  by  separa¬ 
ting  N  into  two  parts: 

N  =  N'  +  N"  , 

where  N'  is  expressed  in  terms  of  the  geoidal  parameters  as  before,  but  where 
N"  should  ideally  represent  the  separation  between  the  qeoid  and  the  measured 
sea  surface. 
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In  this  analysis,  N"  is  expressed  in  terms  of  the  11  tidal  constitu¬ 
ents  considered.  A  physically  meaninqful  model  can  be  conceived  by  allowing 
some  or  all  of  the  tidal  amplitudes  and  phase  angles  to  adjust  within  the 
overall  adjustment  of  satellite  altimetry.  Due  to  a  high  degree  of  uncertain 
ty  in  the  long-period  constituents  (see,  e.g.,  the  "empirical  factor"  mention 
ed  in  [B]),  the  model  for  these  constituents  should  be  linear.  This  will  be 
ensured  by  subjecting  only  the  tidal  amplitudes  to  adjustment,  not  the  phase 
angles.  On  the  other  hand,  since  the  diurnal  and  semidiurnal  tidal  effects 
are  assumed  to  be  known  to  a  better  degree  of  accuracy,  both  these  kinds  of 
parameters  will  enter  the  adjustment. 

Accordingly,  18  tidal  parameters  will  be  added  to  the  adjustment  of 
SEASAT  altimetry:  4  amplitude  parameters  (P.)  due  to  the  long-period  constitu 

J 

ents,  j=A  ,  Mf,  Mm,  SSa;  3  amplitude  parameters  (P.)  and  3  phase-angle  para- 

U  J 

meters  (Q^)  due  to  the  diurnal  constituents,  j*Kj,  Oj,  P^,  and  4  amplitude 

parameters  (P.)  and  4  phase-angle  parameters  (Q.)  due  to  the  semidiurnal  con- 
J  J 

stituents,  j=M,>,  S2,  N^,  K^.  In  considering  that  the  global  adjustment  con¬ 
tains  some  225  terrestrial  parameters  (spherical-harmonic  potential  coeffici¬ 
ents  in  a  14,14  truncated  model)  and  the  regional  adjustment  may  contain  a 
comparable  number  of  parameters  (point-mass  magnitudes),  the  addition  of  18 
parameters  at  either  level  does  not  increase  the  computer  burden  by  a  great 
amount,  both  from  the  run-time  and  storage  requirement  standpoint. 
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4.2  Observation  Equations 


In  agreement  with  Section  4.1,  N",  representing  that  part  of  an 
"instantaneous  geoid  undulation"  sensed  by  the  altimeter  which  can  be  at¬ 
tributed  to  the  tidal  constituents  under  consideration,  can  be  expressed  as 


(4.1) 


N"  =  (hA  +  hMf  +  h^  +  hs$a)  +  (hK  +  hQ  +  hp  )  +  (h„  +  h$  +  hN  +  h„  ) 
0  1112222 


where  the  constituent  heights  h  .  are  the  a  priori  values  (not  equilibrium  in 

J 

general).  The  first  group  identified  by  the  parentheses  represents  the 
long-period  constituents,  the  second  group  represents  the  diurnal  constituents 
and  the  third  group  represents  the  semidiurnal  constituents. 

From  the  above,  the  tidal  part  of  an  altimeter  observation  equation 
is  formed  in  a  familiar  fashion  through  a  linearization  process.  As  indicated 
earlier,  there  will  be  only  one  adjustable  parameter  per  long-period  constitu¬ 
ent,  namely  the  amplitude  parameter  (P^.);  the  diurnal  and  semidiurnal  con¬ 
stituents  will  comprise  two  adjustable  parameters  each  (P.,  Q.).  The  parameter 

I  J 

Pj  represents  a  relative  amplitude  correction  (i.e.,  an  amplitude  correction 

divided  by  the  amplitude  itself)  while  the  parameter  Q.  represents  a  correc- 

1 

tion  in  degrees  to  the  phase  angle.  The  observation  equations  are  developed 
in  terms  of  the  available  a  priori  values  for  h.,  i.e.,  the  linearization  pro- 

J 

ceeds  with  the  initial  values  P.=0,  Q.=0.  When  evaluated  using  the  a  priori 

J  J 

information,  (4.1)  represents  the  tidal  contribution  to  the  constant  term  in 
an  altimeter  observation  equation.  Upon  denoting  the  pertinent  coefficients 
by  p.  and  q.  in  a  self-evident  manner,  the  total  tidal  contribution  to  an 
altimeter  observation  equation  can  be  symbolized  by 


jii  3 


^oV  PMfPMf  +  %PNii  +  pSSaPSSa}  qK1QK1+  P01P01  +  q01°01 

+  PP1PP1+  qp1QP1)  +  (PM2PM2+  qM2QM2+  PS2PS2+  qS2QS2+  PN2PN2+  qN2QN2 
+  p  P  +  q  Q  )  +  N"  ,  (4.2) 


where  the  superscript  "a”  indicates  an  adjusted  value  and  where  N"  is  evaluated 
by  (4.1)  as  stated  above. 

As  is  apparent  from  Section  4.1,  an  original  altimeter  observation 
equation  (without  the  tidal  part)  should  now  be  modified  by  adding  -N"  to  the 
original  constant  term,  augmenting  the  row  of  coefficients  by 

[  “PA  *  -pMf  ’  "  pSSa  ’  'PK  *  ~qK.  ’  "p0  ’  "PP’  "qP.i 

o  111111 

■pM2’  ‘qM2’  ~PS2’  -qS2’  'PN2*  _qN2’  _qK2^’  (4,3) 

and  augmenting  the  column  of  parameters  by 


£  PA0’  PMf ’  PMm’  PSSa’  PK1’  QKj’  P01*  PP1’  °P1; 

pm2>  qm2’  ps2’  qs2’  pn2*  qn2’  pk2’  qk2  ^  * 


(4.4) 


The  formula  (4.2)  can  also  be  written  as 


N,,a  =  l  ha 
j  J 


where  the  adjusted  constituent  heights  are 


The  coefficients  p.  for  the  long-period  constituents  will  be  derived  first 

J 

(the  coefficients  q.  are  immaterial  in  this  case),  followed  by  a  derivation 

*J 

of  the  coefficients  p.,  q.  for  the  diurnal  and  semidiurnal  constituents. 

J  J 

Long- period  constituents.  The  model  equation  for  these  constituents 
corresponds  here  to  the  simplified  model  of  Section  3.1;  in  agreement  with 
the  latter  and  with  [B] ,  Section  3.3,  we  can  write 


h.  =  A.  cosai  , 

J  J  J 

Aj  -  C  Aj  • 
c  =  (1  +  k)e  , 


(4.6a) 

(4.6b) 

(4.6c) 


where 


flj  ■ 


a 


J 

1+k 

e 


equilibrium  amplitude  , 

Greenwich  or  local  argument  (longitude  is 
irrelevant  here)  , 

1.29  (k  is  a  Love  number)  , 

adopted  "empirical  factor"  . 


From  (4.6a)  one  has  at  once: 


dh.  =  dA.  cosa.  =  (A.  cosa. ) (dA./A. )  =  h.P. 

J  J  J  J  J  J  J  JJ 


Accordingly, 


h^  =  h.  +  dh. 

J  J  \J 


Vi +  hj  ■ 


(4.7a) 


and  thus 


Pi  =  hJ 


(4.7b) 
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The  next  step  consists  in  evaluating  h..  Under  idealized  conditions 
(equilibrium  tide  allowing  for  the  earth's  deformation,  see  page  49  of  [B]), 
in  conjunction  with  (4.6c)  one  would  have 

e'  «  1,  c  =  1  +  k  -  1.29  . 

The  value  used  by  Lisitzin  [1974]  --  but  attributed  to  another  source  — 
corresponds  to 

e"  »  3.7  ,  c  -  (1  +  k)  *  3.7  =  4.8  , 

which  was  also  adopted  in  [B] .  However,  the  above  value  for  e"  is  highly 
uncertain  and  could  be  too  high.  To  be  on  the  conservative  side,  we  may 
adopt  an  average  value, 

e  =  |(e'  +  e")  *  2.35  , 

and  thus 


c  *  (1  +  k)  x  2.35  «  3.03  .  (4.8) 

For  the  equilibrium  amplitudes  we  have  from  Section  2.1: 

Ai  =  0.1073  m  (1-3  si n2cf>)  , 

Ao 

A^f  =  0.0131  m  (1  -  3  sin2i}))  , 

=  0.0125  m(l-  3  sin2!}))  , 

A^sa  =  0.0097  m  (1  -  3  sin2!}))  . 

If  these  results  together  with  c  in  (4.8)  are  utilized  in  (4.6  a,b),  one 
obtains 
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^  =  0. 325  m  (1  -  3  si n24>) 

o 

9 

(4.9a) 

w  =  0.040  m  (1  -  3  sin24>) 

C0SaMf  ’ 

(4.9b) 

^  =  0.038  m  (1  -  3  sin2<j>) 

COSotMm  ' 

(4.9c) 

SSa  =  0.029  m  (1-3  sin2<p) 

COSaSSa' 

(4.9d) 

The  arguments  a.  corresponding  to  SEASAT  observational  epoch  appear  in 

J 

(2.42  a-d)  while  a',  appear  in  (2.39  a-c). 

J 

Diurnal  and  semidiurnal  constituents.  According  to  (3.49)  and 
(3.50),  the  model  equation  is 

h-  =  c'£.  +  h  *  (equilibrium  tide).  ,  (4.10) 

J  J  J 

where 

c'  =  0.94  ...  1^  constituent  , 

c‘  =  1.00  ...  the  other  diurnal  and  semidiurnal  constituents; 


£.  =  a.cosa.  +  b.sina.  ...  ocean  tide  ,  (4.10') 

V  J  J  J  J 

a.  =  A.cosiK,  b.  =  A.simJj.  ...  values  computed  from 

J  J  J  J  J  J 

spherical-harmonic  tidal  coefficients  as  in 
(3.52  b,c)  , 

cij  =  Greenwich  argument,  as  in  (2.43a)  -  (2.44d)  with  aj  in 
(2.40a)  -  (2. 41d); 


h  ~  0.62  (a  Love  number)  ; 
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(4.10" ) 


(equilibrium  tide).  =  A’,  cosa.  , 

J  0  3 

A1  =  equilibrium  amplitude,  as  presented  below, 

J 

a.  =  local  argument,  as  in  (2.45  a-g). 

J 


The  A',  for  the  diurnal  constituents  are  gathered  from  Section  2.2 

J 


A'  =  0.125  m  sin  2<j>  , 
K1 

AJ.  =  0.081m  sin  2 <p  , 
U1 


Ap  =  0.047  m  sin  2$ 
K1 


For  the  semidiurnal  constituents  we  have  from  Section  2.3: 


=  0.252  m  cos 2 , 
=  0.113m  cos24>  , 


=  0.049  m  cos2<p  , 


=  0.023  m  cos2<p  . 


Since  (equilibrium  tide),  is  not  adjustable,  one  has 

0 


dhj  *  c' 


From  (4.10')  and  the  relation  that  followed  we  obtain 


£.  =  A.  cosip.  cosa.  +  A.  simp.  sina.  , 

3  3  3  3  3  3  3 


dC-  =  (A.  cosip.  cosa.  +  A.  simp.  sina.)(dA./A.) 

3  3  3  33  3  3  J  J 

+  (l/p)(-A.  simp.  cosa.  +  A.  cosip.  sina. )p  dip.  ; 

J  J  J  J  J  J  J 


(4.11) 
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this  is 


d£.  =  (a.  cosa.  +  b.  sina.)P.  +  (l/p)(-b-  cosa.  +  a.  sina.)0.  ,  (4.12) 

J  J  JJ  J  J  J  J-J  J  J 


where 


pj  5  dW  Qj  E  p  d*j E  d*j  ■ 

with  p= 57.295780  transforming  the  correction  to  the  phase  angle  from  radians 


to  degrees.  Accordingly,  (4.10),  (4.11),  (4.12)  yield 


(4.13a) 


h?  =  c'(a.  cosa.  +  b.  sina.JP.  +  (c'/p)(a.  sina.  -  b.  cosa.)0.  +  h.  , 

J  J  J  J  J  J  JvJJJvJJ 


and  thus 


c' (a.  cosa-  +  b.  sina.)  , 

J  J  J  J 

(c'/p)(a.  sina.  -  b.  cosa.)  , 

J  J  vj  J 

p.  +  h  *  (equilibrium  tide).  . 

J 


(4.13b) 

(4.13c) 

(4.13d) 


We  may  insert  a  remark  that  if  needed,  a  simple  iterative  process 
in  A.  and  iK  could  be  conceived  along  the  following  lines.  From  the  initial 

J 

values  a.  and  b.,  implied  by  the  spherical-harmonic  tidal  coefficients,  one 

J 

computes 


*i  •  <aj  *  bF 


<l»j  =  tg(bj /a. )  . 


The  adjusted  values  can  then  serve  as  new  starting  values  denoted  by  primes: 


a;  -  Aj  ( 1  +  p  j )  , 

*j  '  *j  +  Vp  • 

from  which  the  new  starting  values  and  bi  are  computed  as 

J  J 

aj  =  Aj  cos*j  ■ 

b'.  =  A1,  si m|)'.  . 

3  3  3 


(4.14a) 

(4.14b) 

(4.15a) 

(4.15b) 


If  the  primes  are  omitted  the  adjustment  process  may  start  at  this  point,  as 
if  such  initial  values  were  directly  available. 


4.3  Solution 


The  observation  equations  consist  of  rows  A.. ,  where  i  runs  through 
all  the  observations,  and  of  the  correspondino  constant  terms  .  The 
matrix  of  observation  equations  (A)  consists  of  all  the  rows  A^  and  the  vector 
of  constant  terms  (L)  consists  of  all  the  terms  .  The  altimeter  observa¬ 
tions  are  weighted  independently;  all  of  the  A.  and  l.  are  normalized 
through  the  division  by  the  appropriate  sigma  (square  root  of  the  variance). 
Thus  the  normal  equations  read 

(ATA)X  +  ATL  =  0  ,  (4.16) 

where  X  is  the  vector  of  all  the  parameters.  The  parameters  involved  in  the 
short-arc  adjustment  of  satellite  altimetry  are  divided  into  three  groups: 

1)  corrections  to  the  spherical-harmonic  potential  coefficients,  2)  tidal 
parameters  and  3)  corrections  to  the  state  vector  parameters.  The  first 
group  is  envisioned  as  consisting  of  225  parameters  corresponding  to  a  (14, 

14)  truncated  spherical -harmonic  model.  The  second  group,  absent  in  previous 
adjustments  of  satellite  altimetry  at  AFGL,  comprises  the  18  parameters 
described  earlier.  The  first  two  groups  of  243  "terrestrial  parameters"  are 
assigned  permanent  storage  in  the  computer  core.  The  third  group  consists 
of  six  state  vector  parameters  per  (short)  orbital  arc;  these  parameters  as 
well  as  the  appropriate  portions  of  normal  equations  are  assigned  reusable 
storage  which  is  one  of  the  main  features  of  the  short-arc  algorithm  described 
in  detail  in  several  AFGL  reports.  With  the  aid  of  this  algorithm,  the  state 
vector  parameters  are  eliminated  from  the  normal  equations  and  are  solved  for 
later  (after  the  solution  of  the  terrestrial  parameters),  arc  by  arc. 
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The  row  A.,  in  this  study  is  composed  of  the  "regular"  entries  as¬ 
sociated  with  the  original  adjustment  of  satellite  altimetry  (without  any 
tidal  considerations),  augmented  by  the  18  elements  shown  in  (4.3).  The 
value  L.  is  obtained  by  adding  the  value  -N"  for  the  observation  point  i  to 
the  "regular"  value  of  the  constant  term,  where 

-N"  =  -l  h.  .  (4.17) 

j  J 

This  can  be  expressed  as 

Li  =  L:  +  LV  ,  (4.18) 

where  L.j  represents  the  "regular"  value  of  the  constant  term  and  LV  (computed 
through  equation  4.17)  is  the  contribution  due  to  the  considered  tidal  effects. 
The  vector  ATL  in  (4.16)  is  computed  as 

ATL  =  I  aT  Li  .  (4.19) 

i 

In  the  short-arc  algorithn,  this  operation  is  performed  for  all  the  observations 
on  one  arc,  after  which  the  elimination  of  the  pertinent  six  state  vector 
parameters  takes  place.  Subsequently,  the  same  procedure  takes  place  for  the 
next  arc,  and  so  on. 

At  this  juncture,  a  useful  optional  feature  can  be  incorporated  in 
the  adjustment  algorithm  by  simple  means.  One  may  ask,  what  would  the  vector 
in  (4.19)  become  if  all  the  LV  happened  to  be  zero?  Clearly,  the  answer  is 

ATL-  =  T  aT  Ll  .  (4.20) 

h  1  1 
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B" 


w 


One  should  note  that  no  dimensions  have  changed,  in  other  words,  we  have 
not  reverted  to  the  "regular"  case  without  the  tidal  parameters.  However, 
by  virtue  of  the  extra  effort  undertaken  in  (4.20)  the  computer  program 
is  able  to  solve  an  adjustment  problem  not  only  the  "new"  way  with  the  tidal 
parameters,  but  also  the  "regular"  way. 

The  key  to  this  feature  is  the  option  to  compute  a  "shadow"  vector 
to  aTl.  .  In  particular,  corresponding  to  each  aTl..  there  would  be  a  "shadow" 

T 

vector  A^L‘.  treated  in  exactly  the  same  manner  as  the  original  vector,  i.e., 
one  would  proceed  to  the  summation  (see  4.19  and  4.20),  to  the  elimination 
of  the  state  vector  parameters,  etc.  The  added  tasks  would  in  no  way  change 
the  operations  applied  to  the  matrix  of  normal  equations.  The  added  computer 
burden  in  terms  both  of  the  run-time  and  storage  requirements  would  thus  be 
minimal.  In  order  to  obtain  the  solution  corresponding  to  the  absence  of 
tidal  parameters,  the  normal  equations  would  be  resolved  using  the  finai 
"shadow"  vector  in  conjunction  with  very  large  weights  (i.e.,  very  small  a 
priori  sigmas)  attributed  to  the  tidal  parameters.  Since  this  "shadow"  vector 
corresponds  to  the  zero  constant  terms  L!j  and  since  the  tidal  parameters  would 
be  held  essentially  fixed  at  their  initial  values  which  are  zero  by  construc¬ 
tion,  any  tidal  consideration  would  be  virtually  eliminated  from  the  overall 
adjustment.  It  is  to  be  noted  that  in  a  global  adjustment  of  satellite  alti¬ 
metry,  the  formation  of  observation  equations  and  the  formation  of  normal 
equations  together  with  the  elimination  of  the  state  vector  parameters  require 
one  or  two  orders  of  magnitude  greater  run-time  than  the  actual  solution  of 
normal  equations;  once  again,  the  computer  burden  needed  in  performing  the 
added  task  would  be  minimal.  The  fruit  of  this  undertaking  consists  in  the 
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possibility  of  comparing  the  results  of  a  "new"  adjustment  (in  which  the  tidal 
parameters  are  attributed  realistic  weights)  with  the  results  of  a  "regular" 
adjustment  without  tidal  parameters.  With  the  possible  exception  of  the  C^g 
coefficient  which,  if  weakly  weighted,  could  be  influenced  by  the  exclusion 
of  the  permanent  tidal  term,  the  adjusted  spherical -harmonic  potential  coef¬ 
ficients  should  compare  well  between  the  two  adjustments,  especially  if  the 
altimeter  data  cover  the  oceans  by  a  grid  of  passes  distributed  over  a  suf¬ 
ficiently  long  time  interval.  In  the  presence  of  tidal  adjustment  the  r.m.s. 
residual  is  expected  to  be  smaller  —  by  an  amount  which  is  of  interest  in 
itself  --  than  the  corresponding  r.m.s.  residual  in  the  absence  of  tidal  ad¬ 
justment.  In  both  cases  it  is  assumed  that  the  (weighted)  state  vector  para¬ 
meters  are  uncontaminated  by  tidal  effects;  clearly,  this  assumption  holds 
better  with  the  tidal  adjustment  than  without  it. 

In  the  closing  paragraph  of  Section  4.2  the  possibility  of  using  an 
iterative  scheme  for  tidal  parameters  was  described.  However,  at  the  present 
such  an  approach  is  not  envisioned  for  applications  in  real  data  reductions. 
First,  in  its  original  form  the  iterative  process  could  be  employed  rigorously 
only  if  the  tidal  parameters  were  not  weighted  (or  else  changes  would  have  to 
be  introduced  in  the  constant  terms  so  that  the  weights  would  be  associated 
with  the  initial  values  throughout  the  iterative  process,  not  with  the  updated 
values).  Second,  even  if  a  number  of  intermediate  results  were  saved  on  tapes, 
the  computer  run-time  would  nevertheless  increase  significantly  since  in  each 
iteration,  all  the  steps  associated  with  the  tidal  parameters  and  with  the 
elimination  of  the  state  vector  parameters  (not  to  mention  the  solution  of 
normal  equations)  would  have  to  be  repeated.  Finally,  since  the  tidal  model 
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is  linear  in  tidal  amplitudes  and  nonlinear  in  phase  angles,  an  iterative 

process  would  be  considered  only  if  large  corrections  to  the  phase  angles 

were  suspected.  However,  there  is  little  likelihood  that  this  will  occur 

on  a  global  scale  if  reasonably  good  starting  values  of  spherical-harmonic 

tidal  coefficients  are  available;  the  likelihood  would  be  greater  for  a 

regional  adjustment.  As  a  precautionary  measure,  one  can  ensure  that  these 

corrections  are  within  acceptable  limits  by  attributing  relatively  small 

sigmas  (e.g.,  5°)  to  the  parameters.  On  the  other  hand,  corrections  to 

the  tidal  amplitudes  may  be  large.  If  desired,  or  if  little  a  priori  information 

exists,  the  P.  parameters  can  be  weighted  loosely  without  any  ill  effects,  at 
0 

least  in  theory.  This  property  is  especially  useful  with  regard  to  the  long- 
period  constituents  whose  amplitude  adjustment  may  shed  light  on  the  "empiri¬ 
cal  factor"  e  discussed  earlier. 

A  few  distinct  situations  associated  with  the  weighting  of  tidal 
parameters  can  be  discerned  at  once: 

1)  Large  a  priori  sigmas  (i.e.,  small  weights)  attributed  to  both 
kinds  of  parameters  (P.,  Q.).  In  this  case,  the  amplitudes 
and  the  phase  angles  are  essentially  free  to  adjust.  However, 
as  cautioned  above,  the  linearization  of  the  tidal  model 
could  be  compromised  by  large  corrections  to  the  phase  angles. 

A  narrower  class  of  adjustment  problems  could  be  considered 
where  the  sigmas  attributed  to  Qj  are  no  longer  large. 

2)  Small  a  priori  sigmas  for  Pj,  Q..  The  amplitudes  and  the  phase 
angles  are  essentially  constant  in  the  adjustment.  This  amounts 
to  correcting  the  altimeter  observations  for  the  (presumably) 
known  tidal  effects  and  foregoing  the  tidal  adjustment. 


3)  Small  a  priori  sigmas  for  P^,  CL  in  conjunction  with  zero 
constituent  heights.  This  corresponds  to  zero  L^'  values 
and  to  the  "shadow"  vector  approach  discussed  previously, 
resulting  in  a  "regular"  adjustment  with  no  tidal  effects 
considered. 

The  above  three  groups  are  special.  A  usual  tidal  adjustment  cor¬ 
responds  to  "realistic"  sigmas  which  could  be  envisioned  as  follows: 

Long-period  constituents  . . .  crp  -  1. 

*3 

This  sigma  is  relatively  large  (a  one-sigma  correction  would  imply  a  magni¬ 
tude  of  the  amplitude  itself),  due  to  the  uncertainty  associated  with  the 
"empirical  factor".  Sigmas  of  1.5  or  2  could  also  be  adopted. 

Dirunal  and  semidiurnal  constituents  . ..  op  -  0.5,  aQ  -  10  . 

The  behavior  of  these  constituents  modeled  through  the  spherical-harmonic 
tidal  coefficients  (as  computed  beforehand  with  the  aid  of  Laplace  Tidal 
Equations)  is  assumed  to  be  reasonably  well  known.  One  could  also  consider 
aQ  =  5  (in  0  )  which,  in  view  of  an  earlier  discussion,  might  offer  some 
comfort  with  regard  to  the  nonlinearity  problem. 
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4.4  Practical  Notes 

Tidal  spherical -harmonic  expansion.  As  explained  in  the  closing 
part  of  Section  3.2,  a.  and  b.  for  the  considered  diurnal  and  semidiurnal 

J 

constituents  are  obtained  through  a  spherical -harmonic  expansion  with  the 
tidal  coefficients  taken  from  [Estes,  1980].  These  coefficients  are  com¬ 
plete  through  the  degree  and  order  (12,12).  One  thus  has  for  the  j-th 
constituent: 

12  n 

a.  =  l  J  (a.  cos  mX  +  b.  sin  mX)P  (sin<f>)  ,  (4.21a) 

J  n=o  ro=o  Jnm  Jnm  n 

12  n 

b.  =  l  l  (c.  cos  mX  +  d.  sin  mX)P  ( si nd>)  ;  (4.21b) 

J  n=o  m=o  Jnm  ^nm  nm 

the  ocean  tide  for  this  constituent  is  then  computed  by  (4.10’).  In  these 
equations,  <J>  and  X  are  the  known  geocentric  latitude  and  longitude,  respective¬ 
ly,  of  the  observation  point  associated  with  an  event  on  a  given  satellite  arc. 
When  computing  the  "regular"  geoid  undulation  in  a  global  (spherical-harmonic) 
altimetric  adjustment  the  following  formula,  adapted  from  equation  (2.15)  of 
[Blaha,  1979],  is  used: 

14  n 

N  =  r*  l  (a/r')n  7  (AC  cos  mX  +  ASnmsin  mX)P  (sincj>)  ,  (4.22) 

o  t  nm  nm  nm  r  y 

n=2  m=0 

where  <J>  and  X  have  the  same  interpretation  as  in  (4.21  a,b)  and  where  the 

other  symbols,  explained  in  Section  2.1  of  the  same  reference,  are  of  little 

consequence  in  this  discussion.  The  similarity  in  form  between  (4.21  a,b) 

and  (4.22)  makes  the  computation  of  a.  and  b.  a  trivial  task;  the  sines  and 

J  J 
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cosines  of  the  multiple  angles  mA  as  well  as  the  associated  Legendre  functions 
Pnm(sin<}>)  have  all  been  computed  (by  recursive  formulas)  for  the  evaluation 
of  (4.22).  This  makes  the  spherical -harmonic  formulation  of  tidal  effects 
particularly  attractive  in  conjunction  with  a  global  adjustment  of  satellite 
altimetry. 

Tidal  arguments.  The  state  vector  parameters  on  each  arc,  stored  on 
a  magnetic  tape,  are  preceded  by  the  arc's  ID.  This  ID  gives  the  epoch  time 
in  terms  of  the  day  number  (in  1978),  hours  and  minutes  (UT  time).  For  all 
the  individual  events  on  the  arc,  the  intervals  AT  =  T  -  T  are 

event  epoch 

given  in  seconds  and  fractions  thereof.  The  most  economical  approach  to 
evaluating  the  needed  Greenwich  arguments  is  to  evaluate  them  for  Tgpoch  and 
to  add  a  correction  due  to  AT  by  multiplying  this  interval  by  the  rate  of 
change  in  a'.  (Vsecond),  gathered  from  Section  2.4.  Thus  we  have 

J 

(event)  =  aMf  (epoch)  +  AT  *  0.00030501°, 

aMm  (event)  =  (epoch)  +  AT  x  0.00015122°, 

aSSa^event^  =  aSSa^epoc^  +  AT  *  0.00002282°; 

av  (event)  =  (epoch)  +  AT  »  0.00417807°, 

l  K1 

an  (event)  =  an  (epoch)  +  AT  »  0.00387307°, 

U1  U1 

aD  (event)  =  aD  (epoch)  +  AT  *  0.00415526°; 
y\  *1 


(event)  =  (epoch)  +  AT  *  0.00805114°  , 

ou  (event)  =  ou  (epoch)  +  AT  *  0.00833333°  , 

(event)  =  (epoch)  +  AT  «  0.00789992°  , 

(event)  =  a.,  (epoch)  +  AT  »  0.00835615°  . 

i>2  •'2 

The  local  arguments  follow,  as  usual,  upon  adding  \  to  the  Greenwich  arguments 
for  the  diurnal  constituents  and  2x  for  the  semidiurnal  constituents  (nothing 
is  added  in  case  of  the  long-period  constituents). 

Regional  tidal  adjustment.  Based  on  the  residuals  from  a  global 
adjustment  of  satellite  altimetry,  a  regional  adjustment  in  terms  of  point 
masses  can  be  performed.  This  approach,  described  e.g.  in  the  AFGL  reports 
[Blaha,  1977,  1979,  1980],  results  in  a  more  detailed  local  geoid.  The  point- 
mass  magnitudes  have  been  the  only  parameters  present  in  the  past  "regular" 
approach  to  this  second-phase  adjustment  where  they  have  served  in  accommoda¬ 
ting,  in  a  least-squares  sense,  the  first-phase  altimeter  residuals.  The 
point-mass  adjustment  could  be  complemented  by  tidal  parameters  in  a  manner 
similar  to  that  explained  earlier.  This  possibility  is  now  briefly  outlined. 
Since  the  only  role  of  the  tidal  parameters  in  such  an  adjustment  is  to 
"overcorrect"  the  values  of  P.  and  Q.  from  the  first  adjustment,  the  same 

O  J 

coefficients  in  the  observation  equations  can  be  used  as  those  computed  pre¬ 
viously.  And  since  these  coefficients  are  stored  on  tape  for  each  observa¬ 
tion,  augmenting  the  observation  equations  in  the  second  adjustment  is  a  trivial 
matter.  These  coefficients  are  applied  with  the  positive  sign  because  in  the 
second  adjustment  they  serve  in  accommodating  the  geoidal  residuals  (minus 


the  altimeter  residuals).  The  column-vector  of  parameters  is  then  augmented 
by  the  18  tidal  parameters  appearing  in  (4.4).  The  constant  terms  of  ob¬ 
servation  equations,  stored  on  tape,  are  simply  minus  the  altimeter  residuals 
from  the  first  adjustment.  The  weights  associated  with  the  tidal  parameters 
are  suggested  to  be  similar  to  those  from  the  first  adjustment,  or  perhaps 
greater  since  a  general  fit  of  the  modeled  sea  surface  to  the  observed  sea 
surface  has  been  accomplished  and  is  not  expected  to  change  drastically  in 
the  second  adjustment  if  the  point-mass  region  is  reasonably  large.  Due  to 
the  simple  treatment  of  the  tidal  part  in  the  point-mass  adjustment,  the  run¬ 
time  requirements  are  virtually  unchanged  from  its  "regular"  counterpart. 
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5.  CONCLUSION 


In  the  recent  past,  SEASAT  altimeter  data  together  with  spherical- 
harmonic  potential  coefficients  (supplied  by  the  GEM  10  model)  and  sets  of 
state  vector  parameters  (supplied  by  the  NSWC  precise  ephemeris)  have  been 
adjusted  at  AFGL  through  the  short-arc  algorithm.  In  using  several  low 
degree  and  order  truncations  of  the  spherical-harmonic  (S.H.)  model,  it  has 
been  observed  that  the  empirical  variance  for  geoid  undulations  is  signifi¬ 
cantly  lower  than  the  theoretical  variance.  The  details  of  this  analysis  can 
be  found  in  Appendix  1,  whose  main  findings  are  the  following: 


1)  The  SEASAT  data,  the  ephemeris,  the  S.H.  potential  coef¬ 
ficients  and  the  reference  field  parameters  all  appear  to 
be  of  excellent  quality; 

2)  The  short-arc  algorithm  in  conjunction  with  the  seven- 
minute  criterion  described  in  previous  reports  does  not 
introduce  appreciable  errors;  and 

3)  Within  the  narrow  band  of  frequencies  considered  (cor¬ 
responding  to  the  S.H.  truncations  between  8,8  and 
16,16),  the  adjusted  degree  variances  for  geoid  undula¬ 
tions  reach  only  about  56%  of  their  theoretical  values. 


During  past  efforts  at  AFGL,  the  altimeter  residuals  obtained 
from  the  global  altimeter  adjustment  (with  the  S.H.  potential  coefficients 
and  the  state  vector  components  as  parameters)  have  served  in  regional 
modeling  of  short-wavelength  geoidal  features,  as  well  as  in  studying  geo¬ 
physical  phenomena  such  as  ocean  bottom  topography.  In  this  adjustment 
the  geoid  has  been  assumed  to  coincide  with  the  ocean  surface  as  sensored 
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by  the  altimeter  and  the  sea-surface  effects  have  been  ignored.  However, 
due  to  increasing  geophysical  interest  in  a  realistic  representation  of  the 
open  ocean  tide,  the  latest  development  of  the  short-arc  satellite  alti¬ 
metry  model  allows  for  the  inclusion  of  the  most  important  tidal  constitu¬ 
ents.  In  particular,  an  adjustment  algorithm  has  been  designed  in  which 
four  long-period  constituents,  three  diurnal  constituents  and  four  semi¬ 
diurnal  constituents  may  be  subject  to  adjustment  within  the  overall  adjust¬ 
ment  of  SEASAT  altimetric  observations.  Except  for  the  long-period  con¬ 
stituents  where  the  phase  angle  is  considered  fixed,  both  tidal  amplitude 
and  phase  angle  are  adjustable.  This  development  is  contained  in  the  body 
of  the  present  report. 

One  of  the  building  blocks  of  the  tidal  adjustment  just  mentioned 
is  the  exploitation  of  the  S.H.  expansion  of  the  constituent  height.  In 
particular,  the  S.H.  tidal  coefficients  for  the  diurnal  constituents  Kj,  Oj, 
Pj  and  for  the  semidiurnal  constituents  M^,  S2,  Ng,  can  be  adopted  from 
[Estes,  1980],  or  from  another  reference  on  the  subject,  and  can  be  treated 
with  advantage  by  essentially  the  same  algorithm  which  serves  in  the  compu¬ 
tation  of  geoid  undulations  from  the  S.H.  potential  coefficients.  One  has 
to  keep  in  mind,  however,  that  very  serious  errors  would  occur  if  the  tidal 
coefficients  were  truncated  (the  above  reference  supplies  a  set  of  12,12 
coefficients  for  each  of  the  seven  constituents).  This  property  stems  from 
the  method  used  in  the  determination  of  these  coefficients,  in  which  the 
land  areas  were  ignored  (the  solution  was  not  held  to  zero  over  land).  Thus, 
as  stated  on  page  80  of  [Estes,  1980],  "completely  spurious  values  are  ob¬ 
tained  if  the  expansion  is  evaluated  over  land  areas."  This  has  been  con¬ 
firmed  by  evaluating  a  few  constituent  heights.  For  example,  although  for 
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the  middle  of  the  U.S.  none  of  the  constituent  heights  exhibit  an  unreason¬ 
able  magnitude,  the  M2  tide  for  the  middle  of  Asia  (latitude  45°,  longitude 
90°)  at  a  given  time  epoch  (day  200  in  1978,  UT = 2  hours)  is  computed  as  7.1m. 

In  addition  to  exercising  proper  caution  when  dealing  with  the 
tidal  coefficients,  one  should  view  the  individual  tidal  constituents  in 
relation  to  the  specific  orbital  characteristics  of  the  satellite  used  in 
data  acquisition.  In  particular,  the  constituents  which  cannot  be  separated 
from  others  by  the  altimeter  data  should  be  adequately  modeled  and  enforced 
a  priori.  If  left  to  the  adjustment,  such  constituents  would  result  in 
aliasing  tidal  constituents  of  different  frequencies.  For  example,  if  a 
satellite  performing  altimeter  measurements  were  sun- synchronous  (this  is 
not  the  case  with  SEASAT),  the  constituents  and  S2  would  be  aliased  into 
the  zero  frequency,  i.e.,  this  particular  satellite,  sampling  theoretically 
only  the  high  solar  tide,  would  sense  them  as  a  constant  constituent.  In 
general,  then,  one  is  faced  with  the  necessity  of  attributing  a  higher  weight 
to  the  parameters  (here  the  tidal  amplitude  and  the  phase  angle)  of  each 
tidal  constituent  which  cannot  be  satisfactorily  resolved  from  the  altimetry. 

In  the  present  adjustment  model  of  SEASAT  altimetry  and  the  selected 
tidal  effects,  the  constituents  to  be  treated  with  care  are  Kp  P^  S2  and  K^. 
As  stated  in  [TOPEX,  1981],  the  diurnal  constituents  K1  and  are  aliased  to 
a  six-month  and  a  constant  constituent,  respectively;  the  semidiurnal  con¬ 
stituents  S2  and  l<2  are  similarly  aliased  to  a  six-month  and  a  three-month 
constituent.  Since  the  useful  life  span  of  SEASAT  amounted  to  only  about 
three  months,  most  of  these  constituents  cannot  be  properly  resolved,  as  is 
already  the  case  with  SSa  (the  semiannual  constituent).  To  further  complicate 
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the  matters,  Kj  and  S^,  if  not  partially  or  totally  constrained  by  appropri¬ 
ate  weighting,  would  be  aliased  into  this  SSa.  Accordingly,  it  may  be  neces¬ 
sary  to  lower  the  a  priori  sigmas  associated  with  the  eight  parameters  cor¬ 
responding  to  the  above  four  constituents  (K^,  Pp  and  K^,  two  parameters 
per  constituent),  perhaps  to  %  or  less  of  the  sigmas  associated  with  the 
other  diurnal  and  semidiurnal  constituents. 


APPENDIX  1 


ANALYSIS  OF  LOW  DEGREE  VARIANCES  FOR  GEOID  UNDULATIONS 
USING  SEASAT  ALTIMETER  DATA 

The  global  adjustment  of  SEASAT  altimeter  data  yields,  at  an 
initial  stage,  the  misclosures  in  the  altimeter  observation  equations 
formulated  in  terms  of  spherical-harmonic  (S.H.)  potential  coefficients 
and  state  vector  parameters  (six  per  short  orbital  arc).  These  misclosures 
are  also  called  the  constant  terms  of  observation  equations.  They  reflect 
on  the  errors  in  the  altimeter  measurements  (the  noise  of  the  system),  in 
the  a  priori  state  vector  parameters,  in  the  short-arc  algorithm  (a  model¬ 
ing  error),  and  in  the  input  S.H.  potential  coefficients  and  reference  field 
parameters.  But,  most  of  all,  they  reflect  on  the  geoidal  detail  ignored  by 
the  adjustment.  Such  detail  depends  on  the  truncation  of  the  S.H.  model 
and  can  be  represented  by  a  variance,  called  here  the  "theoretical  variance", 
obtained  as  a  sum  of  "theoretical  degree  variances"  for  geoid  undulations, 
where  the  summation  extends  over  all  the  neglected  degrees.  Thus,  in  a  model 
truncated  at  degree  and  order  (14,14)  the  summation  extends  over  the  degrees 
15  through  perhaps  1,000  (beyond  this  degree  the  contribution  is  negligible 
for  all  practical  purposes).  The  theoretical  degree  variances  are  in  turn 
computed  from  the  covariance  function. 

When  an  actual  adjustment  with  SEASAT  data  was  carried  out,  it  was 
observed  that  the  root  mean  square  (rms)  of  the  altimetric  misclosures  in 
the  (14,14)  truncated  S.H.  model  was  significantly  lower  than  the  "theoretical 
sigma"  (the  square  root  of  the  theoretical  variance).  Yet,  this  rms  contains 
error  contributions  from  the  other  sources  of  information  mentioned  above, 

-74- 


J 


w 


i.e.,  from  the  ephemeris,  the  S.H.  potential  coefficients  together  with  the 
reference  field  parameters,  and  the  SEASAT  altimetry.  The  low  rms  value 
serves  as  an  indication  of  excellent  quality  of  these  sources  and,  in  addi¬ 
tion,  as  an  indication  of  sufficient  accuracy  in  the  short-arc  algorithm 
applied  in  conjunction  with  the  seven-minute  arc  criterion  discussed  in  pre¬ 
vious  reports.  This  is  supported  by  the  fact  that  the  average  misclosure  is 
nearly  zero  (-0.1m).  The  rms  value  suggests  that  the  theoretical  formula 
for  the  covariance  function  may  be  too  conservative,  at  least  insofar  as  the 
geoid  undulations  for  relatively  low  degree  and  order  truncations  are  con¬ 
cerned. 

In  analyzing  this  phenomenon  the  variance  due  to  the  ephemeris,  in 
particular,  due  to  the  radial  component  of  the  state  vector  parameters,  has 
been  subtracted  from  the  mean  square  of  the  misclosures  yielding  the 
"empirical  variance".  This  implies  that  the  other  variances  affecting  the  rms 
(the  variances  due  to  the  input  S.H.  coefficients  and  the  reference  field 
parameters,  to  the  altimeter  noise,  and  to  the  short-arc  algorithm)  have  been 
assumed  to  be  negligible;  otherwise  the  already  low  empirical  variance  would 
have  been  even  lower.  The  adopted  procedure  for  computing  this  variance  is 
symbolically  expressed  as 

a2  -  rms2(misclos>'es)  -  az(ephemeris)  ,  (Al.l) 

where  "A"  represents  "empirical".  An  initial  step  in  the  misclosure  analysis 
was  undertaken  in  Appendix  1  of  [Blaha,  1981]  where  the  sources  of  information 
beyond  the  SEASAT  altimeter  were  identified  (GEM  10  for  the  S.H.  potential 
coefficients,  NSWC  precise  ephemeris  for  the  state  vector  parameters,  and 
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IUGG/IAG,  1975,  for  the  reference  gravity  field  parameters).  The  S.H.  model 
considered  therein  was  limited  to  the  (14,14)  subset.  However,  the  global 
rms  misclosures  have  been  since  computed  also  in  conjunction  with  the  (8,8), 
(10,10),  (12,12),  and  (16,16)  truncated  models.  In  order  to  assess  the 
usefulness  of  these  results  for  further  analysis,  the  topic  of  the  theoretical 
variance  will  be  addressed  in  more  detail. 


If  an  adjustment  process  is  carried  out  in  terms  of  an  (n,n)  S.H. 
model,  the  theoretical  variance  for  geoid  undulations  characterizing  the 
higher  frequency  content  ignored  by  the  adjustment  is  obtained  as  a  summation 
of  theoretical  degree  variances  from  degree  n+1  onward.  This  variance  is  ex¬ 
pressed  as 


a 


2 


n+l,<*> 


I 


k=n+l 


» 


where  c°  is  the  k-th  degree  theoretical  variance  for  geoid  undulations.  The 
symbol  »  is  replaced  in  practice  by  a  suitable  large  number  such  as  1,000. 
According  to  the  formulas  in  [Tscherning  and  Rapp,  1974]  the  theoretical  k-th 
degree  variance  for  geoid  undulations  (cr£)  is  computed  from  the  theoretical 
k-th  degree  variance  for  gravity  anomalies  (a£)  which  is  in  turn  obtained  from 
the  covariance  function,  as  follows: 

o>  =  R2  52/[G2(k-l)2]  , 

52  =  sk+2A(k-l)/[(k-2)(k+B)]  , 


where  G  is  the  average  value  of  gravity  (979.8  gal),  R  is  the  earth's  mean 
radius  (6,371  km),  and 


-76- 


s  =  0.999617  , 

A  =  425.28  mgal 2  , 

B  =  24,  an  integer  . 

These  formulas  have  been  adapted  from  equations  (10),  (25A)  and  from  Table  7 
of  the  same  reference. 

Upon  carrying  out  the  operations  indicated  above,  the  theoretical 
k-th  degree  variance  for  geoid  undulations  reads 

a*  =  0.999617k+2  17 ,98 lm2/ [( k- 1 ) ( k-2 ) ( k+24)]  .  (A1.2) 

Upon  summing  up  these  variances  from  the  chosen  degree  on,  one  obtains 
the  results  presented  in  the  second  column  of  Table  1.  The  "n+1"  in  the 
first  column  implies  the  (n,n)  truncated  S.H.  model  and  indicates  that  the 
geoidal  variances  (both  theoretical  and  empirical)  are  due  to  the  neglected 
degrees,  from  n+1  onward.  The  empirical  variances  are  featured  in  the  third 
column.  They  are  computed  according  to  (Al.l),  where  the  rms  misclosures  are 
obtained  from  the  altimeter  observation  equations  using  the  (n,n)  truncated 
model ,  and  where 

a(ephemeris)  =  1.6  m, 

which  is  the  a  priori  standard  error  in  the  radial  component  of  the  state  vector 
parameters.  It  is  apparent  from  Table  1  that  the  empirical  variance  for  the 
low  degree  truncations  considered  is  significantly  lower  than  its  theoretical 
counterpart.  For  the  truncation  (8,8)  the  ratio  of  these  two  variances  is 
0.4884,  for  the  truncation  (10,10)  it  recedes  to  0.4004  and,  eventually,  for 
the  truncation  (16,16)  it  rises  to  0.4917.  Since  a  similar  pattern  can  be 
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assumed  also  for  the  intermediate  truncations,  one  can  conclude  that  for 
any  truncation  between  (8,8)  and  (16,16)  this  ratio  maintains  itself  under 
0.5.  However,  it  would  most  likely  climb  above  0.5  already  for  the  truncation 
(17,17)  due  to  the  rising  uncertainty  in  the  S.H.  potential  coefficients  with 
the  increasing  degrees.  In  accordance  with  the  outlined  procedure,  this 
effect  should  ideally  be  eliminated  from  the  empirical  variance  (the  vari¬ 
ance  due  to  the  input  S.H.  coefficients  was  previously  assumed  to  be  negligible). 


1 


i 


n+1 

an+l,“ 

an+l ,°° 

57.758 

28.21 

4U.758 

16.32 

12+1 

30.562 

12.88 

14+1 

23.886 

10.13 

16+1 

19.243 

9.46 

Table  1  | 

Theoretical  variances  (in  mz)  and  empirical  variances  (in  m2)  ! 

for  geoid  undulations  corresponding  to  selected  degree  and  order  truncations  (n,n) 


The  marginal  usefulness  of  the  S.H.  model  beyond  (14,14)  for  the 
study  of  the  empirical  variances  is  corroborated  by  the  statement  on  page  19 
of  [Khan,  1981]:  "At  this  stage  it  is  generally  believed  that  the  long  wave¬ 
length  components  of  the  gravity  field  are  wel 1 -determined  to  n=14".  Figure  2 
of  this  reference,  showing  the  approximate  level  of  accuracy  as  a  function  of 
frequency,  quantifies  this  statement  for  the  GEM  10B  coefficients  through  the 
degree  n=22  by  showing  that  for  n=14  the  "accuracy  determination"  is  54%  while 
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for  n=16  it  is  only  about  33%  and  it  continues  to  diminish  as  n  increases. 
When  faced  with  these  facts  one  could  choose  to  disregard  the  empirical 
variance  for  n=16  and  perhaps  even  for  n=14,  and  use  only  the  first  three 
empirical  variances  in  Table  1.  However,  in  considering  that  only  the  dif¬ 
ferences  between  the  empirical  variances  will  eventually  be  used  as  obser¬ 
vables  in  evaluating  the  formula  (A1.2),  it  becomes  clear  that  such  an  ap¬ 
proach  would  drastically  reduce  the  already  small  number  of  available  ob¬ 
servations. 

Another  avenue  one  could  take  in  analyzing  (A1.2)  would  be  to  use 
all  of  the  empirical  variances  listed  in  Table  1  and  to  correct  them  for  the 
effect  of  uncertainty  in  the  input  S.H.  coefficients.  However,  several  dif¬ 
ficulties  would  complicate  this  approach: 

a)  The  accuracy  evaluation  in  Figure  2  of  [Khan,  1981]  concerns 
the  GEM  10B  coefficients  and  not  the  GEM  10  coefficients 
which  have  served  in  the  computation  of  empirical  variances 
from  SEASAT  altimetry; 

b)  For  a  given  degree  k,  the  "accuracy  determination",  y^%,  is 
not  specific  enough  to  be  unequivocally  applicable  to  the 
problem  at  hand.  In  particular,  it  is  not  clear  whether 
the  effect  of  the  uncertainty  in  the  k-th  degree  S.H.  coef¬ 
ficients  is  exactly  { 1-y^ )%  of  the  theoretical  degree 
variance  a£,  or  some  ether  percentage;  in  other  words,  it 
is  not  clear  what  probability  statement  could  be  associated 
with  the  uncertainty  (l-yk)o£  ; 

c)  Even  if  (l-yk)ak  were  known  to  be  the  theoretical  effect  of 

a£,  rather 

than  (l-yk)ak,  should  be  subtracted  from  each  uncorrected 


this  uncertainty  one  would  suspect  that  (l-yk)t 
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empirical  degree  variance;  here  t  is  some  reducing  factor  due 
to  the  smaller  magnitudes  of  the  empirical  variances  as  com¬ 
pared  to  their  theoretical  counterparts;  and 

d)  If  an  adjustment  model  for  degree  variances  should  be  repre¬ 
sented  by  the  formula  (A1.2)  or  its  equivalent,  and  if  the 
observables  should  be  represented  by  the  corrected  empirical 
degree  variances  or  their  combinations,  one  would  still  be 
faced  with  the  task  of  attributing  proper  weights  to  the  ob¬ 
servations  in  order  to  carry  out  a  reasonably  rigorous  least- 
squares  adjustment  process;  for  the  sake  of  simplicity  (but 
without  any  further  justification),  one  might  wish  to  adopt 
a  unit  weight  matrix  for  this  purpose,  but  the  variance- 
covariance  matrices  after  adjustment  would  be  difficult  to 
interpret  in  a  meaningful  manner. 

In  the  approach  adopted  herein  the  empirical  degree  variances  are 
not  corrected  in  any  way,  but  the  uncertainty  in  the  S.H.  coefficients  serves 
in  the  establishment  of  a  weighting  scheme  in  the  least-squares  adjustment 
with  (A1.2)  as  the  basic  model.  The  latter  is  rewritten  as 

a*  =  sk+2  C/[(k-l)(k-2)(k+B)]  ,  (A1.3) 

where  the  parameters  are  s,  C  and  B.  The  initial  values  of  these  parameters 
were  presented  earlier  as 


s 

=  0.999617  , 

(A1.4a) 

c 

=  17,981  m2  , 

(A1.4b) 

B 

=  24. 

(A1.4c) 
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Since  various  approximations  applied  to  a  weighting  scheme  are  much  less 
dangerous  than  the  same  kind  of  approximations  carried  into  the  actual  ob¬ 
servations,  the  difficulties  encountered  in  the  preceding  paragraph  are 
now  unimportant.  By  the  same  token,  the  association  of  the  measure  of  un¬ 
certainty,  (l-yk)a£  ,  with  the  standard  error  (sigma)  of  the  quantity  o£ 
would  be  acceptable  if  this  standard  error  were  indeed  necessary.  However, 
if  the  actual  observations  are  the  empirical  degree  variances  (ajp  which, 
for  the  degrees  involved,  are  systematically  lower  than  the  theoretical 
degree  variances,  one  can  take  as  the  a  priori  standard  error: 


sigma  (ajp  =  (l-yk)t 


(A1.5) 


where  the  reducing  factor  t  is  computed  as  the  ratio  between  the  sum  of  the 
pertinent  empirical  degree  variances  (here  between  the  degrees  9  and  16)  and 
the  sum  of  the  corresponding  theoretical  degree  variances, 

t  =  °9, 16  1  a9, 16  ' 

The  degree  variances  have  been  limited  to  n=9  through  n=16  for 
practical  reasons.  As  stated  in  [Blaha,  1981],  the  errors  introduced  in 
the  short-arc  algorithm  due  to  the  degree  and  order  truncation  below  (8,8) 
cannot  be  tolerated  in  SEASAT  altimetry.  In  fact,  even  the  truncations  (8,8) 
and  (9,9)  introduce  detectable  (but  tolerable)  errors  in  the  satellite  posi¬ 
tions,  but  the  set  (10,10)  is  already  completely  satisfactory.  It  should 
be  noted  that  these  considerations  are  valid  only  if  the  arc's  duration  does 
not  surpass  7  minutes  and  if  the  state  vector  parameters  are  given  at,  or 
very  nearly  at, the  mid-arc.  For  these  reasons  the  empirical  variances  in 
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Table  1  are  listed  starting  with  n+1-8+1.  On  the  other  hand,  the  last 
value  listed  corresponds  to  n+1  =  16+1,  due  to  the  large  uncertainty  in  the 
S.H.  coefficients  discussed  earlier.  Clearly,  with  the  observations  avail¬ 
able  only  for  the  degrees  indicated,  the  results  of  the  present  analysis  are 
applicable  only  to  the  low-frequency  range  of  the  S.H.  model,  between  n=9 
and  n=16.  Since  the  desired  combinations  of  degree  variances  are 


16 


J9,16 


£ 


k=9 


J9 ,16 


16 


k=9 


2  2 

Oq  -  a;-,  , 

9,°°  17, »  * 


^  2  ^z 

a „  -  ai-,  , 

9,“  17,“  ’ 


Table  1  yields 


t  =  (28.21m2-  9.46m2)/ (57.758m2  -  19.243  m2)  *  0.4868  .  (A1.6) 

The  variance  of  the  quantity  a £  is  computed  from  (Al.5)  as 

var(a2)  =  [(l-yk)t  a£]2  .  (A1.7) 


For  the  weighting  purposes,  the  empirical  variances  for  any  degrees  are  con¬ 
sidered  mutually  independent.  Since  the  empirical  variances  will  be  considered 
in  combinations  (in  fact,  they  will  be  added  pair-wise,  each  pair  forming  one 
observation),  one  has 

n2 

var(a2  n  )  =  l  var(S2)  .  (A1.8) 

1’  2  k=n-^ 


The  diagonal  weight  matrix  will  thus  be  composed  of  the  reciprocal  values  in 
(A1.8).  Although  these  variances  are  only  very  approximate  --  and  will  be 
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further  simplified  below  --  their  magnitudes  are  realistic  and  allow  meaning¬ 
ful  estimates  of  the  variance-covariance  matrices  after  adjustment. 

Based  on  the  foregoing  discussion,  one  can  compute  the  desired 
weights  upon  adopting  the  accuracy  determination  factor  (denoted  here  as  yk%) 
from  Figure  2  of  [Khan,  1981].  This  factor  is  listed,  in  the  form  (l-y^X, 
in  the  third  column  of  Table  2  for  k  between  9  and  16.  In  the  cases  k=9  and 
k=10  the  values  (l-y^  are  augmented  by  7%  and  3%,  respectively,  due  to  the 
errors  introduced  in  the  short-arc  algorithm  via  the  truncation  (8,8).  Al¬ 
though  these  percentages  would  be  far  too  high  for  the  present  configuration 
of  SEASAT  short  arcs,  a  conservative  approach  is  warranted  because  the  pertin¬ 
ent  rms  misclosures  were  computed  before  the  completion  of  the  new  short-arc 
preprocessor  which  ensures  that  the  epoch  is  almost  exactly  at  mid-arc.  The 
fourth  column  of  Table  2  gives  the  variances  of  according  to  (A1.7),  with 
t  taken  from  (A1.6).  Since  each  of  the  observations  will  consist  of  two 
empirical  variances  added  together,  the  observational  variances  are  formed 
according  to  (Al. 8),  where  n2=ni+^>anc*  are  listed,  for  the  degrees  9  and  10, 

11  and  12,  13  and  14,  15  and  16,  in  the  fifth  column  of  Table  2.  It  is  ap¬ 
parent  that  these  variances  (in  m4)  are  fairly  close  to  unity.  Since  only 
approximate  weights  are  needed  and  since  several  approximations  in  their  for¬ 
mation  have  already  been  introduced,  little  damage  will  be  done  if  the  weight 
matrix  is  adopted  as  the  unit  matrix. 
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k 

<'-V* 

var(a2) 

var(S£bn2) 

9 

10 

9.689 

7.311 

5%+7%=12% 

12%+3%=15% 

0.3203 

0.2850 

0.61 

fl 

5.680 

4.516 

19% 

27% 

0.2760 

0.3523 

0.63 

B 

3.661 

3.015 

36% 

46% 

0.4116 

0.4558 

0.87 

B 

2.517 

2.126 

57% 

67% 

0.4878 

0.4808 

0.97 

Table  2 

Variances  (in  m1*)  of  the  empirical  degree  variances 
and  of  their  pair-wise  combinations 


As  has  already  been  indicated,  each  observation  equation  corresponds 
to  two  degree  variances.  The  constant  terms  in  these  equations  are  formed 
as  the  sum  of  two  theoretical  degree  variances  minus  the  sum  of  two  corres¬ 
ponding  empirical  degree  variances.  The  pair-wise  sums  of  the  theoretical 
degree  variances  are  apparent  from  the  second  column  of  Table  2.  They  can 
also  be  obtained  from  the  second  column  of  Table  1  as  the  differences  between 
the  appropriate  theoretical  variances;  for  example,  erg  jq  is  obtained  (in  m2) 
from  Table  2  as  9.689  +  7.311  =  17.000  or  from  Table  1  as  57.758  -  40.758  = 
17.000.  The  corresponding  combinations  of  the  empirical  degree  variances  are 
apparent  from  the  third  column  of  Table  1,  for  example,  a g  1Q  =  28.21  -  16.32  = 
11.89  (in  m2).  Table  3  lists  the  pertinent  sums  of  theoretical  and  empirical 
degree  variances  (columns  two  and  three,  respectively) ,  and  of  their  differences 
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L 


under  the  heading  L 
to  the  differences 


nl’n2*  tota^s  ^or  co^umns  two  and  three  correspond 
"top-bottom"  in  the  same  columns  of  Table  1. 


1  9  ^2 

a2  n 
nl’n2 

BIB 

Ln1,n2 

9,10 

17.000 

11.89 

5.110 

11,12 

10.196 

3.44 

6.756 

13,14 

6.676 

2.75 

3.926 

15,16 

4.643 

0.67 

3.973 

I 

38.515 

18.75 

19.765 

Table  3 

Constant  terms  Lp^p^  (in  m2)  in  the  observation  equations 


The  matrix  of  observation  equations  in  the  case  of  individual  degree 
variances  (see  A1.3)  would  contain  the  partial  derivatives  of  with  respect 
to  s,  C  and  B: 

9o2/9s  =  a2(k+2)/s  , 

9cr£/9C  =  cr£/C  , 

9a2 /9B  =  -a2/(k+B)  . 


-85- 


The  observation  equations  for  the  sums  of  degree  variances  read 


=  oJ(k+2)/s]ds  +  n  ^/C]dC-[I  cr£/(k+B)]dB +  L  . 

1  2  k^n-^  k=n^  1  2 


In  order  to  make  the  numbers  in  the  matrix  of  observation  equation  (A)  com¬ 
parable,  one  may  replace  ds,  dC  and  dB  by  ds',  dC‘  and  dB '  in  such  a  way 
that 


ds  = 

(1/100)  ds*  , 

(A1.9a) 

dC  = 

(C/10)  dC*  , 

(A1.9b) 

dB  = 

-4dB'  ; 

(A1.9c) 

the  observation  equations  then  become 

n2  n2 
Vn  n  =  [  l  0.01  ajJ(k+2)/s]ds'  +  [  l  0.1cr*]dC' 

1’  2  k=n2  k^j^ 

n2 

+  [  l  4a//(k+8)]dB '  +  L  n  .  (A1.10) 

k=n  *  nl’n2 


The  elements  of  the  matrix  A  formed  from  (A1.10)  are  computed  in  a 
straightforward  fashion  using  the  values  in  (A1.4  a-c)  and  from  the  second 
column  of  Table  2  (the  additions  in  A1.10  are  pair-wise  as  before): 


1.9438 

1.7000 

2.0345 

1.3712 

1.0196 

1.1509 

1.0320 

0.6676 

0.7132 

0.8109 

0.4643 

0.4708 
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(Al.ll) 


The  observation  equations  are  written  in  the  matrix  form  as 

V  =  AX  +  L  , 

where  V  is  the  vector  of  (four)  residuals,  X  is  the  vector  of  (three)  para¬ 
meters,  ds',  dC'  and  dB',  and  L  is  the  vector  of  (four)  constant  terms  as 
presented  in  Table  3. 

Since  the  weight  matrix  has  been  stipulated  to  be 
P  =  I  , 

where  the  units  are  m"4,  the  least-squares  solution  is  given  by 

X  =  -N-1  AT  L  , 

N  =  ATA  , 

and  the  variance-covariance  matrix  for  the  parameters  X  is 

cx  =  N_1* 

The  original  parameters  (denoted  by  Y)  are  computed  as  in  (A1.9  a-c).  From 
the  general  relation 

Y  =  GX  , 

Cy  =  G  Cy  G^, 

it  follows  for  the  specific  linear  functions,  such  as 
*1  =  Clxl  * 

y2  -  c2x2  ’  ••• 
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that 


IG1 '  %  • 


Jy2  |c2‘  °x2  •  ••• 


and 


Pv  v  =  si9n  (cic-)  P  , 

y1.y2  1  t-  xi»x2 


yielding  the  sigmas  and  the  correlation  coefficients  (p)  for  the  original 
parameters  ds,  dC  and  dB. 

As  may  have  been  anticipated  already  from  (Al.ll)  showing  that  the 
values  in  all  three  columns  of  A  behave  in  a  quite  similar  manner,  the  matrix 
of  normal  equations  (N)  is  severely  ill-conditioned.  In  fact,  both  X  and  Cx 
contain  exceedingly  large  and  unrealistic  values.  The  matrix  of  correlation 
coefficients,  formed  from  N_1  and  modified  so  that  the  signs  correspond  to 
the  original  parameters,  reads 


(A1.12) 


1.  -0.999079  -0.998350 

p  =  J  symm.  1.  +0.999892 

1.  J 

Although  the  results  of  this  adjustment  are  meaningless,  the  rms  residual  has 
been  computed  for  comparison  purposes  as 


rmSy  =  0.668m2. 

This  value  would  be  larger  if  dB  were  rounded  to  the  nearest  integer. 
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Clearly,  all  three  parameters  cannot  be  determined  in  any  meaning¬ 
ful  way  from  the  observations  that  encompass  only  the  degrees  n=9  to  n=16. 

One  could  attempt  to  adjust  only  two  parameters  while  holding  the  third 
fixed  (the  corresponding  column  in  A  would  be  disregarded).  However,  since 
any  two  of  the  three  columns  in  (Al.ll)  are  quite  comparable,  ill-conditioning 
is  again  a  problem.  For  example,  in  terms  of  the  original  parameters  s,  C  one 
obtains 


ds  =  -  0.09618  ,  a$  =  0.02730  , 

dC  =  +13.880  m2  ,  aQ  =  6,225  m2. 

Both  the  corrections  and  the  sigmas  are  quite  large.  The  post-adjustment  cor¬ 
relation  coefficient,  although  lower  than  any  one  in  (A1.12),  is  still  very 
large,  namely 

PsC  =  -0.990871  . 

the  residual  fit  is  necessarily  worse  than  previously  (with  three  parameters), 
represented  by 

rmsv  =  0.874mz. 

The  adjustments  of  the  other  two  pairs  of  parameters  is  equally  unimpressive; 
in  fact,  the  correlation  coefficient  for  dC  and  dB  reaches  +0.998928. 

Another  possibility  of  adjusting  two  parameters  is  to  restrict  the 
model  to  only  one  adjustable  parameter  (s,  C  or  B)  and  to  add  a  linear  para¬ 
meter,  D.  The  following  results  are  obtained  for  C  and  D: 
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dC  =  -2,058  m2  ,  oc  =  1,914  m2  , 

D  =  -1.920  m2  ,  aD  =  0.570  m2  ; 

pc0  =  -0.899  ; 

rmSy  =  1.017  m2  . 

Comparable  results  would  be  obtained  also  if  dC  were  replaced  by  ds  or  dB 
(if  the  latter  were  rounded  to  the  nearest  integer  the  rms  residual  would 
necessarily  increase).  However,  it  is  not  the  purpose  of  the  present 
analysis  to  change  the  structure  of  the  degree  variance  model.  Accordingly, 
this  example  has  served  merely  as  an  illustration  of  a  less  ill-conditioned 
system  (on  the  other  hand,  the  rms  residual  has  increased  somewhat). 

It  has  been  already  noticed  that  due  to  the  increased  uncertainty 
in  the  S.H.  potential  coefficients,  the  truncation  (16,16)  may  be  of  only 
marginal  value.  For  this  reason  another  series  of  adjustments  has  been  car¬ 
ried  out,  but  without  the  last  observation  equation.  All  the  results,  in¬ 
cluding  the  goodness  of  fit  represented  here  by  the  rms  residual,  have  suf¬ 
fered  when  compared  to  their  counterparts  with  four  (rather  than  three)  ob¬ 
servation  equations.  The  worsening  is  most  noticeable  in  the  correlation 
coefficients.  A  certain  amount  of  similarity  has  been  observed  only  when 
adjusting  one  parameter  alone.  The  results  of  the  one-parameter  adjustment 
will  now  be  presented  for  the  original  four  observation  equations.  With  s 
as  the  parameter,  one  obtains 

ds  =  -0.035862  ,  og  =  0.003681  ; 

rmSy  =  1.417  m2  . 

By  comparison,  the  adjustment  of  B  yields 

dB  =  +14.751  ,  oB  =  1.607; 

rmSy  =  2. 163  m2  . 
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This  residual  fit  would  be  even  worse  if  dB  were  rounded  to  15.  One  notices 
that  the  adjusted  parameter  Ba=  B  +  dB  would  be  39  which  is  greatly,  perhaps 
unrealistically,  removed  from  the  theoretical  value  B  =  24. 

The  most  useful  in  the  series  of  adjustments  carried  out  in  the 
course  of  this  analysis  is  probably  the  adjustment  of  the  single  parameter 
C.  Upon  obtaining  the  adjusted  parameter  Ca  =  C  +  dC,  one  can  construct  the 
ratio  Ca/C  as  a  simple  measure  of  judging  the  decrease  in  degree  variances, 
theoretical  versus  empirical,  within  the  range  n=9  and  n=16.  The  results 
obtained  with  this  model  are 

dC  =  -7,849  m2  ,  ac  =  839  m2  ; 

Ca  =  10,132  m2  ,  (Al. 13a) 

Ca/C  =  0.5635  .  ( Al. 13b) 

The  result  (A1.13b)  indicates  that  within  the  narrow  band  of  the  S.H. 
frequencies  treated  herein  the  adjusted  degree  variances  would  reach  only 
about  56%  of  their  theoretical  value  (the  adjusted  sigmas  would  accordingly 
reach  about  75%  of  their  theoretical  value).  The  rms  residual  is  quite  high 
in  this  case,  namely 

rms v  =  1.966  m2  . 

This  suggests  that  due  to  the  unaccounted  for  errors  in  the  S.H.  potential 
coefficients  and  perhaps  to  other  causes,  the  empirical  degree  variances  do 
not  behave  exactly  according  to  the  model,  the  correcting  factor  Ca/C  not¬ 
withstanding.  But  an  important  conclusion  that  can  be  reached  is  that  in 
spite  of  these  errors  the  empirical  variances  for  the  considered  (low) 
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frequencies  are  substantially  smaller  than  their  theoretical  counterparts, 
and  that  the  ratio  of  these  two  kinds  of  variances  is  approximately  0.56. 
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APPENDIX  2 

EXPANSIONS  IN  SURFACE  SPHERICAL- HARMON ICS 

Most  of  the  formulas  used  in  this  appendix  can  be  found,  in  one 
form  or  another,  in  the  standard  literature  and  will  thus  be  stated  without 
proof  or  reference.  In  the  derivations,  use  will  be  made  of  the  spherical 
coordinates  0  (colatitude)  and  X  (longitude). 

A  function  on  the  sphere,  f(8,X),  can  be  expanded  in  surface 
spheri cal -harmoni cs : 

f(9,X)  =  l  f  (0,X)  ,  (A2.1) 

n=0 

where  the  surface  spherical-harmonic  of  n-th  degree  is 

fn(6-x>  -  S„  Kn,  +  bnm  Snmle-X>l  -  I*2'2’ 

m=0 

a's  and  b's  being  the  coefficients,  and  R's  and  S's  being  defined  as 

Rno(9’A)  =  Pno(cos9)  3  pn(cos0)  .  Sn0{9’A)  =  0  » 

Rnn](0,X)  =  pnm^cos0)  cos  mX  »  m>0  , 

Snm(0,X)  =  Pnm(cos9)  sin  mX  ,  m>0  , 


where  pnm(cos0)  are  the  associated  Legendre  functions  in  the  argument  cos©. 
Equation  (A2.1)  can  thus  be  rewritten  in  a  familiar  form  as 


(A2.3) 


f(9.A)  =  l  I  [a 

n=0  m=0 


R  (0,X)  +  b  S 
nm  nm  nm  nm 


(6,X)J 


oo  n 

=  Y  J  (a„m  cos  mX  +  b„„  sin  mX)P„m(cos0)  . 
~  nm  nm  nm 

n=0  m=0 
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If  the  spherical  angle  separates  two  points  with  the  spherical 
coordinates  (0,X)  and  (6,X),  the  following  "addition  formula"  applies: 

(A2.4) 

Pn(cos<j,)  =  Pn(cos0)  P  (cose)  +  2  l  [(n-m)!/(n  +  m)!][R  (e,X)  R  m(0,X) 

m=l  nm  nm 

*  s«»<0-x>  • 

where 

cos  ip  =  cose  cose  +  sine  sine  cos  (I-  X)  . 

Next,  the  orthogonality  relations  are  presented: 

-«./<*♦  1)  •  (A2.5.) 

(A2.5b) 

ff[Rnm{0^)J2dfi  =  J5[Snm(0,X)]2dQ  =  [2ir/(2n  +  1)]  [(n  +  m)  !/(n  -  m) !]  , 

m>0  , 

where  dft  is  the  solid  angle  element  and  n  represents  the  surface  of  a  unit 
sphere.  It  is  important  to  note  that  the  integrals  of  all  the  other  possible 
products  not  covered  by  (A2.5  a,b)  vanish;  one  type  of  such  products  is 

*k£<0*X>-  ^  ^  n  and/or  £  i  m,  where  the  R's  could  be  replaced  by 

the  S's,  and  another  type  is  Rnm(0,X)  S^fe.X)  for  any  n,  m,  k,  t. 

The  relation  (A2.3)  leads  to  the  statement 

I(6,x)  -  tt  f(e.D  P„(cosrt  d n-ssj'  RM(e,X)  *  bM  s M(e.x)] 

«  P  n ( c  os  ^ )  dn  , 

associated,  as  indicated,  with  the  (fixed)  point  (0,X).  Upon  expanding  Pn(cos\|>) 
according  to  (A2.4)  and  utilizing  the  orthogonality  relations,  one  obtains 
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I(9’A)  =  Rno(9’X)  anot47T/(2n  +  1)]  +  2  j?  [(n  -  m)  !/(n  +  m)  !l 

m=l 


anm[2Tr^2n+l)]  [(n+m)  !/(n-m) !]  +  s  (e,A)  b 

nm'  nm 

*  E27r/(2n+l)]  [(n+m)  !/(n-m)  !j  } 

This  last  result  is  written  at  once  as 


!«•.»>■  [«./(*♦«]  jo  &nm  Rnm(e,x)  .  biw  s^e.x,]  , 


Which,  with  the  aid  of  (A2.2),  yields  the  desired  relation: 
fff(Q,X)  Pn(cos^)  dfi  =  [4tt / ( 2 n  +  1)] f  (e,A)  . 


(A2.6) 
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